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Preface 


The purpose of this thesis has been to present a self- 
seneained text explaining the fundamentals of the high altitude 
electromagnetic pulse in a manner which would provide a basic 
understanding of the phenomenon for the first year graduate 
student of Engineering. 1 would like to emphasize the fact. 
that it has not been the purpose of this thesis to provide 
"state of the art" threat levels for the fields of the 
electromagnetic pulse. Raters the calculations concerning 
these levels which are included in this report are meant 
only to show order of magnitude agreement of the model 
developed here with accepted theories. Although some of 
the results of derivations in the report involve advanced 
electromagnetic theory, I have attempted to develop these 
derivations in a manner that requires only a familiarity 
with integral and vector te Maxwell's equations, and 
a few basic precepts of modern physics. 

I would like és express my Sincere appreciation to 
Major Carl T. Case for the many hours of guidance and advice 
which he devoted to this effort. I would also like to thank 
Dr. Charles J. Bridgman for his encouragement and interest 
in the project. Finally, I would like to acknowledge my 
wife's patience and support, without which I could not have 


completed this thesis. 
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Abstract 


This thesis constitutes a microscopic approach to the 
siadeeousgneere pulse (EMP) caused by a high altitude nuclear 
weapon burst. This approach involves the detailed calculations 
for the time-dependent electromagnetic fields of a single 
relativistic, Compton electron created in the upper atmosphere 
by the gamma radiation from the weapon detonation. The effect 
of ail Such electrons created is superimposed to obtain the 
net effect or the EMP measured by an observer. The basic 
theory developed ignores the conductivity which builds up 
in the surrounding air. This is shown to be a relatively 
good assumption for weapons with gamma radiation yields less 
than a kiloton. A brief section is included to indicate the 


derivations necessary for the EMP in a conducting media. 
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THEORY OF THE HIGH ALTITUDE 


ELECTROMAGNETIC PULSE 
I. Introduction 


With the advent of nuclear weapons several previously 
unsuspected phenomena emerged as possible system damage 
mechanisms. One such phenomenon is the interaction of the 
radiation from a nuclear weapon with the environment in 
such a manner as to cause the motion of electric charges. 
The motion of these charges in turn causes a pulse of Sioa 
tric and magnetic fields which can be very intense and can 
be radiated for great Airstanees: The nature of these fields 
varies according to the altitude at which the weapon is 
detonated. This thesis will concern only those fields 
associated with weapon detonations which occur at exoatmos- 
pheric altitudes. With the development of weapon systems 
involving highly sensitive electronics, the electric and 
magnetic fields of the electromagnetic pulse (EMP) asso- 
ciated with high altitude bursts constitute a phenomenon 


which should be fully understood by defense planners. 


Definition 

The EMP can be defined as a long range nuclear weapon 
effect wherein a high intensity, padioetcaqueney electro- 
magnetic pulse is generated (Ref 9, chap X). There are 
several other similar nuclear effects which should not be 
confused with EMP. These include the Argus effect, the 


blackout phenomenon, and transient effect on electronics 
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| 
(TREE) . The Argus effect is the phenomenon whereby the 
electrons produced by a nuclear burst are trapped in the 
earth's magnetic field thereby creating a potential radia- 
on hazard for personnel and electronic equipment (Ref 2, 
pg 1). The blackout phenomenon relates to the ionization 
produced by a nuclear detonation in or above the atmosphere 
and the subsequent radio frequency interference (BLACKOUT) 
feecieing from the attenuating, reflecting, cluttering and 
scattering of radar and radio signals (Ref 2, pg 1). TREE 
refers to the direct effects of radiation on electronic ‘ 


systems and components (see also Ref 14, pp 1-3). 


History 

Although it has been said that Enrico Fermi predicted 
the EMP before the Trinity explosion on 16 July 1945 (Ref 14, 
pg 4), the high frequency component of this signal associated 
with high altitude weapon @atonetione did not become a seri-_ 
ously studied problem until the early 1960's. There are 
several reasons for the recent interest in this problem. One 
of these is the fact that the electronics in weapons systems 
have become more sensitive to extraneous electrical fields. 
Furthermore, it has been shown that the high altitude EMP 
can serve as a "signature" for a weapon detonation, thereby 
allowing for greater control of the ban on atmospheric test-: 
ing. Also, there has been a greater understanding of the 
‘impact of the effect on systems previously thought hardened 
by defense planners. 


Some of the first serious studies of the problem of the 
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high altitude EMP began in 1961 with theoretical studies by 
W. J. Karzas and R. Latter ( Ref 12).. These studies pro- 
vided theoretical models and mechanisms for describing what 
was then the unconfirmed high frequency component of the high 
altitude EMP. The existence of this component was confirmed 
in the 1962 Pacific tests in which "a component of the 

"radio flash' which depends on the magnetic bearing of the 
nuclear detonation from the observation point" was observed 
(Ref 20, pg 3). From the early 1960's until the SrAceEe 

the EMP, and more specifically the high altitude, high 


frequency EMP, has received a great deal of theoretical 


and experimental study. 


High Altitude EMP Source 

The EMP which results from a high altitude weapons burst 
is generated by electrons which are created in the upper at- 
mosphere. These electrons are created by energetic photons 
(y-rays) from the weapon which undergo collisions with tne 
molecules of the air, ejecting atomic electrons. These 
scattered electrons are very energetic and radiate by means 
which will be described in a later chapter. 

It should be noted that electromagnetic signals can 
arise from a nuclear burst from sources other than that 
mentioned above. An example of one such signal is the eae 
signal which is associated with assymetries in the weapon 
casing. Another involves the electromagnetic fields gen- 
erated within a system by the interaction of the weapon 


radiation with the materials of the system itself (System- 
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Generated EMP or SGEMP). A third example concerns the 
electrostatic fields in the near vicinity of the burst 
cused by charge separation due to ionization (Close-In 
EMP). The relative threats of these other effects are small 


as compared to the high altitude EMP. 


Microscopic Versus Macroscopic 

When calculating field strengths associated with the 
“high altitude EMP, the general approach is "macroscopic" 
in nature. That is, one solves Maxwell's equations for 
electromagnetism using expressions for the "current" of 
electrons created in the upper atmosphere. Calculations 
are Made based on the "current density" or the amount of 
charge passing through a unit area in a unit time interval. 
If one can find an expression for the current density, then, 
in principle, one can compute the electric and magnetic 
fields associated with this current. 

Although the approach described above has been used 
with much success to obtain actual values for field strengths 
associated with the EMP, it is not an enlightening approach 
for one who is interested in the physical processes which are 
involved in the generation of the EMP. There is a consid- 
erable amount of physics “buried” in this macroscopic approach 
to the EMP. For this reason, it is the purpose of this thesis 
to develop a microscopic model of the high altitude EMP. 
Rather than calculating the EMP field strengths based on a 
current density, the fields for a single, accelerated electron 


have been calculated in great detail. The effect of all the 
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electrons created are then superimposed to obtain the over- 
all or net field strength of the generated EMP. It is hoped 
that this microscopic approach will enable the serious stu- 
dent of the EMP to have a more basic and satisfying under- 


standing of the high altitude EMP. 


Overview 

The first chronological step in the generation of the 
high altitude EMP is the production of gamma rays resulting 
from a nuclear burst. These gamma rays+travel through the 
upper atmosphere depositing energy through Compton collisions 
thereby creating energetic electrons (See Figure 1). In 
Chapter II, which deals with the nature of the weapon gamma 
radiation and its transport through the atmosphere, it is 


shown that the number density of Compton electrons, N at 


pri’ 
a @istance R from a burst occurring at an altitude of HOB 


is given by 


exp{u Ze 108/20 (1-e®/%o) } 


(II-17) N__.(R,t) = A(R) 
Boe 41R? 


t 
x {gr fecer- 2) att} 
Y oO 


where A is the Compton macroscopic cross baewieas Uo is the 
total macroscopic cross section at sea level, Zo is the 
scale height, xy is the total gamma yield, EY is the average 
gamma ray energy, f(t) is the time dependence of the yield, 
and an exponential atmosphere has been assumed. 

The Compton electrons are turned by the earth's 


magnetic field which causes them to radiate. It is shown 
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in Chapter III that, for a given observer, the radiated: 
fields which he measures result from only those electrons: 
moving within a small cross-sectional area centered about 
his line-of-sight with the weapon burst (see Figure 1). 
Further, it is shown that, for the periods when its radiation 
is significant, the electron hee essentially a constant speed. 
The power radiated by such an electron is primarily in the | 
forward direction (see Figures 15 and 16) as is indicated 
by the following formula (developed in Chapter III) for the 
power radiated per solid angle (8): : 
aP 


; a ae 1 2 
(III-40) ag- = Gee Ga) (I-68 cos jp)? ia - ¥2(1-8 cos y)? 


sin*y cos? n 
where C, E, @ are the speed of light, free space per- 
mittivity, and electron charge respectively, 8 if the 
electron's relative speed, 8 is its relative acceleration, 
and w and n are angles defined in Figure 14. If A(R) 
represents the cross-sectional area which encloses electrons 
contributing to the EMP measured by an observer located at a 


distance X from the burst, then the time dependent fields 


which he observes are given by: 


-- (KER) 


Cc 
A(R) A(R) F(R) i att E(2") Bi (ter'") 
oO 


x 
(III-46) E(X,t) = g* far ~“Ah 
oO 


Y 
where F(R) is the R dependence given by the exponential term 
divided by 4n7R? in (II-17) above and E,(t) is the general 
time dependence of the field strengths radiated by a single 
electron and measured by the observer. This equation is 
developed in Chapter III. 

Determining the size of A(R) is one.of the key cal- 


6 
“4 
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culations in this approach. Chapter ‘IV illustrates how to 
obtain this variable.and also gives some results of the EMP 
which an observer would measure. These calculations were 
performed assuming that the conductivity of the air was 
zero. The range of validity for this assumption is checked 
in Chapter V and it is indicated how one might proceed 

when air conductivity is significant. The equations above 
and the developments which follow are relativistically 


correct. 


GNE/PH/74-7 


WEAPON BURST 


; ELECTRONS 
CREATED 
CYCKCICIGRCIC 
“RARTH'S SAG verte’ 
FIELD-OUT OF / 
PAPER rs 


' ELECTRONS TURNED BY 
\BARTE'S wAGNETIC FIELD 
1 \ : 


p=2 
RADIATION PATTERN 
FOR SINGLE ELECTRON: 


eer ne 


SVIRGIN GAMMAS 


Figure i Generation “of the High avectude EMP. 
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II. Review of Nuclear Weapon Physics and Effects 


Since the EMP depends upon the radiation emitted from 
a nuclear weapon, it is the purpose of this chapter, to pre- 
sent a brief discussion of the physics involved in the creation 
and transport of the radiation produced from a nuclear weapon 
jekouseion. 

The chapter will encompass three main areas relating to 
weapon radiation. The first area will involve the origin 
and nature of the radiation that is produced within the - 
initial weapon burst. Secondly, the transport of this 
radiation through the surrounding atmosphere will be con- 
sidered. Finally, the interaction of the radiation with 
the atmosphere to create energetic electrons (Compton effect) 
will be investigated. Key derivations are included in the 
Appendices sna nee referenced throughout the main body of 


the text. 


With regard to the high altitude EMP the important 
portion of the weapon output is the eadiacies in the form 
of gamma rays and x-rays. For weapon detonations occurring 
below 100 or so kilometers, the Seintey source of the EMP 
involves only the electrons created by the gamma rays. One’ 
reason for this is that at these altitudes the range of the 
electrons created by the x-rays is too short to cause an 
appreciable current. Furthermore, for detonations at 
these altitudes, the X-rays are absorbed in a relatively 


small region about the burst point. This region has been 
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made conducting by ionization from the earlier .emission of 

gamma rays and hence the weak radiation from the electrons 

created by the x-rays is effectively damped out (Ref ll, 

pp 40-5). Since the basic physics for the creation of the 

EMP is essentially the same for x-rays and gamma rays, only 
the EMP generated by gamma rays will be considered in this 

thesis. 

There are two main characteristics of the gamma radi- 
ee output of the weapon which must be considered with 
respect to the EMP. First it is adeedenns to know the 
energy of the gamma rays and second, one must know the time 
dependence of the radiation yield of the weapon. The energy 
of the gamma rays produced in a typical theoretical fission 
device range from under 1.0 MeV up to around 5.0 MeV, with 


an average energy of approximately 1.5 Mev. 


For simplicity, 
in any calculations which follow it shall be assumed that the 
spectrum of gamma rays emitted from a weapon is monoenergetic 
such that all gamma rays have an energy of 1.5 MeV. The time 
dependence of the radiation yield of a hypothetical weapon is 
characterized by a sharp rise followed by a decay both of 
which are generally described as exponential in nature. The 
physical reasons for the time dependence of the radiation 


output is beyond the scope of this study, however, the in- 


terested reader is referred to The Physics of Nuclear 


linformation concerning the average energy was obtained 
in private correspondence with Dr. C. J. Bridgman, Physics 
Dept., AFIT. 


10 
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Explosives which is a set of class notes by Dr. C. Jd. 
Bridgman (Ref 1). 

.There have been several attempts to analytically 
describe the time dependence of weapon yield. One such 
attempt which is an empirical model for the time depend~ 
ence presents the following formula. 

(at 8) expfa(t-t.)} 
(II-1) £(t) = N(Bta expt (ars) (E-t,) Fy (Ref 18) 


where a and 86 are weapon parameters, te peeaumacene time 
t at which the time dependence f(t) is m maximum and N is a 
normalization parameter chosen such that 

00 
(II-2) J f(t) at= 2 

ce) 


A graph of this particular time dependence is given in 


Figure 2 using- typical values of o and 8. 


Pee etna ree Semen nares Sinn 


After one has an understanding of the nature fenahey 
and time dependence) of the radiation emitted from a nuclear 
weapon, it is necessary to gain an understanding of how this 
radiation is transported through the surrounding environment. 
That is, if there is a given amount of radiation produced at 
the weapon burst, how much of this radiation will reach some 
point which is a known distance away. 

If the gamma ray pulse emitted from a nuclear deto- 
nation can be described with respect to time by a function 
f(t) as mentioned in the previous section then at any time 
t the number of gammas produced per second is given by (next 
page) . : . 


ll 
i 
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: aN. (t) Y £(t) 
(II-3) Se so ike ea a 


dt EY 
where a is the gamma yield of the explosion (in MeV), EY is 
the average photon energy (in MeV ~ assuming a monoenergetic 
source), and NL (t) is the number of gamma rays. It can be 
shown (see Appendix A) that if the weapon burst is assumed 
to be a point source of gamma's then at a distance R from 
this source the number of unreacted or virgin gammas passing 


through a unit area per unit length of time (the virgin 


intensity - 6) is given by” . 
R 


anitecty 2 pe 
2 eh aca. 
(II-4) ®(R,t) at = gaa 


where u(r) is the total macroscopic cross section (see 
Appendix A) for photon reaction with the environment and c 

is the speed of light. As is evidenced by (II-4), when con- 
sidering high altitude EMP one is concerned only with those 
gammas that arrive at a point R without having been scattered 
or acted upon by the environment in any manner (virgin inten- 
sity - 6). The reason for this is that the high frequency 
EMP associated with high altitude bursts is essentially a 
"line-of-sight" phenomenon. That is, if an observer is below 
the burst, then only those electrons that are traveling 
directly on the "line-of-sight" between the observer and the 
weapon burst will contribute significantly to the EMP that 


the observer measures (see Figure 3). This will be more 


rhe s in the argument on Ny accounts for time delay 


due to transport of y - pulse over distance R. 


13 
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Figure 3. EMP as Line-of-Sight Phenomenon. The fields from 
the electron created at (B) will contribute to the 
EMP measured at (C) whereas the electron created 
‘at (D) will make no contribution. 


fully explained in Chapter III. Although -it is possible that 
a gamma photon which was not originally on the "line-of-sight" 
could be scattered onto the Tipdseeienee thereby coestaae 
an electron which could possibly contribute to the EMP, the 
chances for this are small enough that it will be ignored 

in the rest of this development. 

Alenousn (II-4) provides a formalistic solution to the 
problem of the transport of virgin gamma photons, it is de- 
sirable to have a form of ©(R,t) with which one can obtain 
Mepiad chee eal values. The main problem with this in (II=4) 
is the following integral in the exponential term: 


R 
(II-5) fur) ar 


u | _ 
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In order to obtain a closed solution of this integral, it 

is necessary to make .some approximations concerning the 
nature of the function u(r). First it shall be assumed 

that the macroscopis cross section u(r) varies with r only 
as the density of the media through which the gamma photons 
are being transported varies with r. In air, a good approx- 
imation is to let the denarey vary exponentially with the 


altitude Z above sea level (See Figure 4) such that 


(II-6) — p(Z) = Poe ° ‘ . 


where 0(2) is the density of air at an altitude Z, oe is 
the density of air at sea level, and Zo is a scaling factor 
(if 2 is measured in kilometers then Z%7)- If one lets Uy 
represent the macroscopic cross section at sea level and 


makes the substitution Z = HOB-r then using the approximations 


outlined above u(r) becomes 


(II-7) u (r) a Uo (HOB-r) /Z, 


The integral in (II-5) can now be evaluated using (II-7) such 


that 
R R - (HOB-r) /Z 
(II-8) JS u(r)dr = f Hoe ° ar 
re) fe) 


Upon integrating, the right hand side of the above becomes © 


R ~ (HOB~r)/Z, ~HOB/Z r/Z R 
f ule dr = uZe (e 
re) 


Oo 
oO 0 o 
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Figure 4. Geometry of High Altitude EMP. 


or ; 
Se «ee  -HOB/2Z, R/Z , 
(II-9) Sf wlr)ar = wo2oe- Or Ga “244 (3) 
: (e) 


Substituting (II-9) into (II-4) one has a closed solution 


for the virgin intensity of gamma photons at a distance R 


such that ; 
R -HOB/Z | R/Z, 
GE) 2 RiRe) See Eee, oor 
| - 47R? 
aN 


¥Y given in (II-3) the virgin 


t 


Using the expression for 


3the approximations made in solving this integral were 
obtained through private correspondence with Capt. Terry 
Cc. Chapman, GNE-74M, AFIT. : 


l 
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intensity can be written as 


-HOB/Z R/Z 
- Y exp{u ze (l-e )} 
(II-11) ®(R,t) = = Fo en eae a ae te 
Y 4nR? 


The reader should be aware that it has been tacitly assumed 
that the macroscopic cross section u(r) is independent of 
time. For early time periods following a weapon burst this - 
is a safe assumption. The air chemistry involved in late 
time calculations concerning u(r) are beyond the scope of 


this thesis. 


Creation of Electrons 

There are several ways in which gamma photons can inter- 
act with the molecules of air to produce electrons. These 
processes include the photoelectric effect, pair production, 
and the Compton effect. The probability of any one of these 
events occurring to a given gamma photon depends upon the 
energy of that photon and the material through which it is 
traveling. This fact is illustrated in Vigure 5 on the 
following page which is a graph of the mass attenuation 
ccefficients* corresponding to the various processes as a 
function of the energy of the photon in air. The curve 
corresponding to total attenuation is the overall sum of 
all the individual effects. With respect to Figure 5 one 
can see that for relatively low energies (<0.1 MeV) the 


photoelectric effect is the predominant factor in the total 


4ohe mass attenuation coefficient is the macroscopic 
cross section (see Appendix A) divided by the density of 
the photon environment. 
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Photoelectric 
Absorption 


Total 
Absorption 


% 


Compton ee ~~ 
Absorption %5 be 
Production——t 
roduction 
LP 
0.0001° 
0.01 


Fig. 5. Mass Attenuation Coefficients for 
Photons in Air. The curve marked “total 
absorption" is (“%) = 4) +(e) +p) 
where Oc , 7 and « are the corresponding linear 
coefficients for Compton absorption, photoelectric 
absorption, and pair production; “© is the air 
density. For example at 1.5 Mev the Compton 
absorption is 2.4 x 10°75 meter?/kg of ‘air; at 
0° C and 760 mm Hg pressure, the density of air 
is Pf = 1.293 kg/m3. (Refi14, p,40) 
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attenuation whereas for relatively high photon energies 
(>10 MeV) pair production provides the major contribution. 
However, for the intermediate photon energies (between 0.1 
and 10 MeV) the Compton effect is the predominant inter- 
action for photon attenuation. Therefore, since weapon 
gamma radiation consists of photons with energies between 
1 MeV and 5 MeV the Compton process is the effect which 
shall be considered with respect to the creation of electrons. 
The Compton effect (scatter) is the process wheeeby a 
gamma photon is scattered by an atomic electron such that 
the energy of the gamma photon is decreased and the electron, 
which will be assumed to be initially at rest, is ionized 
and set into motion. For gamma photons of 1.5 MeV the 
atomic electrons, because of their relatively low binding 
energy, may be: considered as free electrons which are at 
rest. The geometry of the situation is fiiperdated in 
Figure 6. ‘If EY is the incident photon energy then one. 
can easily show (see Appendix A) that the energy of the 


scatter photon is given by the relation 


; 
(I-12) B= —— 

rae 
where 6 is the angle between the direction of the incident 
and scattered photon (see Figure 6), mo is the electron 
rest mass, and c is the speed of light in a vacuum. The 


kinetic energy of the scattered electron Te can be easily 


found using the conservation of energy (assuming elastic 
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Figure 6. Geometry of Compton Collision. (h is Planck's 
constant and v is the photons' frequency) 


collision). That is 
2 ae ‘ 
(1 cos 6)E,/m.¢ 


a $. = 2 aT =  aray 2 aR WT PARRA: 
(II-13) Ts. Ey E Y ‘EV ligdeos OE 7c? 


As was mentioned earlier, the only electrons itch are of 
interest with respect to the high altitude EMP are those: 
which are moving on the "line-of-sight" that. is, those 
electrons corresponding to @ = 0. This is not entirely 
true. It will be shown that electrons which Gary slightly 
off the "line-of-sight" do contribute to the EMP. Although 
the majority of electrons scattered by 1.5 MeV gamma photons 
are in the "line-of-sight" direction (Z=0), one can show 
that there is a distribution of directions for the scat~- 
tered elections ranging from $=0° to @=90° (Ref 17, pp 16-9). 


It is possible to account for those electrons scattered off 
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the "line-of-sight" but which still contribute. to the EMP. 
However, this is a secondary effect and will be ignored in 
this’thesis. Therefore it will be assumed that all the 
electrons created by the gamma photons travel in the initial 
direction of the gamma photon (#=0°). To find T. for these 
electrons, one uses (EI-13) with 9=180° since momentum must 
be conserved. Thus 


E 

(11-14) 1 oo 
e-LoS 1 + mc?/2E 

o Y 


. . 


where T is the kinetic energy for the "line-of-sight" 


e-Los 
@lectrons. 
Since it will be necessary later to eupenampeee the 
individual effects of these Compton electrons, it is de- 
Sirable to know how many of these are created at various 
distances from the burst. A quantity which is needed in 
determining this is the reaction rate or the number of 
Compton events which, occur in a unit volume over a unit 
time period. If X} is the macroscopic cross section for 
the Compton interaction and $(r,t) is the virgin intensity 
of gamma photons at a distance r from the burst, then the 


reaction rate for Compton scatter of virgin photons is 


given by the relation 
A(II-15) R, (x,t) = A(r) 6(r,t) 


where Ro is the reaction rate. If one assumes that each 
one of these events creates an energetic electron which 
contributes to the EMP (it's motion is on the "line-of-~ 


sight") then the time rate of change of the number density 
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of these "primary" electrons would be given by the relation 


ON ora (Tet) 


at = Baler’) 
or 
(II-16) dn... (r,t) 
; pri. oe 
at A(r) @(xr,t) 
where N_.. is the primary electron number density. Therefore, 


pri 


at any time t after the detonation and at any distance R from 
the point of detonation, the number density of primary elec- 


trons is 


t 
Nppi (Rt) = A(R) S a(R EN) atT 
or using equation (II-10) and (II-4) 


-HOB/Z R/Z 
exp{y_Z_e P(l-e =) } 
(II-17) N nq (Ret) =A (R) 9 -——____________ 

P ATR? 


t 

Y R 

x {a fs £(t' - =) at} 
EY 6 Cc 


To see physically what is represented by (II-17) one 
should understand how the gamma energy is deposited through- 
out the atmosphere with respect to distance from the burst. 


This can be done by considering (II-17) for late times 


(te). If one assumes a normalized time function then in 
the limit as t approaches » the expression for Nori becomes 
-HOB/Z R/Z 
. yy exp{u,2.¢ (l-e )} 
(II-18) N ry (Retr) = RAD eet See eae ee ee eee 
P Y ATR? 


The function A(R) can be approximated using the exponential 


approximation illustrated in equation (II-7). A graph of 
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Nong (Ret) is given in Figure 7 on the following page where 
Uo? which is the total attenuation cross section, has been 
taken to be equal to the cross section for the Compton inter- 
action.” It should be noted that it is a rather common 
practice in much of the technical literature concerning the 
high altitude EMP to use A for Wy in the above expression 
‘since the Compton effect is a predominant interaction in 

the energy ranges of interest. This is seen to be a valid 
assumption if one refers to Figure 5 which shows that the 
total attenuation coefficient is due almost entirely to 
Compton attenuation for 1.5 MeV gamma photons. With regard 
to Figure 7, it is seviade that most of the primary electrons 
are created in a layer of air between thirty and fifty kilo- 
meters. For this reason, when speaking of the high altitude 
EMP, one often refers to the “absorption layer" which is 

that portion of the atmosphere in which most of the gamma 


energy is absorbed and most of the primary electrons are 


created as described above. 


Summary 

There have been several assumptions and points made 
in this chapter which should be reiterated at this time. 
First, although the gamma photons which emanate from a 
nuclear weapon burst have an energy spectrum ranging from 
around 1.0 MeV to 5.0 MeV, it will be assumed for the re- 


mainder of this study that the gamma spectrum is mono- 


*rhese values were obtained from the American 
Institute of Physics Handbook (Ref. 5, p. 8-93). 
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ahergenie with an energy of 1.5 MeV. With respect to the 
‘creation of primary electrons it has been asserted that 
only the virgin or unscattered gamma photons are pen Cane 
for the nigh frequency EMP. Finally, the attenuation of 
the virgin photons with energy of 1.5 MeV has been assumed 
to be due to the Compton effect and all other processes 


have been ignored. 
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III. Electromagnetic Theory of the High Altitude EMP 


The purpose of this chapter is to describe and S develop 
is adenedaopie model. Difficult and unwieldy derivations ap- 
pear in the Appendices and are referenced in the body of the 
chapter. 

The first section of this chapter is devoted to the gen- 
eral field equations as developed from Maxwell's equations for 
a single electron. This is followed by a section concerning 
the equations of motion of an electron in the upper atmosphere. 
The results of the first two sections are then combined in the 
third section to obtain-the field equations for an electron 
‘moving in a specific manner. The effects of all those elec- 
trons which contribute to the high altitude EMP are then 
summed up in the fourth section. The effects of air conduc- 
tivity are ignored in this chapter but are discussed in 


Chapter V. 


General Field Equations for a Single Electron 


The most basic equations in electromagnetic theory are 


Maxwell's equations. They are presented here without justi- 


fication or derivation as follows® 
<= 3B (r 
Vx E(xr,t) = SS SB et) 
V x B(r,t) = HJ (r,t) tule, ob iret) 
(III-1) eee 
se p(t, t) 
Voe« E(r,t) = age 
Vv 
V - B(x,t) = 0 
6 


The interested reader is referred to the text Electro- 
magnetic Fields and Waves by Lorrain and Corson. 
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where E is the electric field, B is the magnetic field, Pe is 


the charge density, J is the current density, We is perme- 
ability of vacuum, e, is the permittivity of vacuum, r is 
the position at which the fields are measured, and t is the 
time at which they are measured. 

If one defines the quantities A(r,t) and V(r,t) such 


that’ 


_ dA(x,t) 
at 


E(r,t) = -vv(r,t) 


(III-2) ae ae : 
B(r,t) = V x A(r,t) ‘ 
and then manipulates the equations (III-1) using the above 


substitutions, the following relations are obtained (see 


Appendix B): i 
a a Tere p(e,t) 
(III-3) V2v(r,t) + aE {VeA(r,t)} = - —— 
ey 
_A an 4 _ o7A(r t) Ses er 
(III-4) V*A(r,t) Uyey aes We, o(r,t) 
. t 
ok oe ; ovi(r,t) 
+ V{VeA(r,t)} + Viue, oer a } 


Since a vector field is not uniquely defined unless one 
specifies both its curl and divergence and since only the 
curl of A has been specified in the above, one can "choose" 
the following relationship concerning the divergence of A 


av(r,t) 


(III-5) VeA(r,t) = - Hoey, aE 


This procedure is known as picking the guage of A and the 
particular gauge chosen above is called the Lorentz gauge. 


Upon substituting (III-5) into (III-3) and (III-4) one 


"RK is often called the vector potential and V the scalar 
potential. 
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obtains the following wave Agaations: 


a ; 2 =o Pp. (re) 
vev(z,t) - 2 Svteet) _ - c 
- ce? ot? Vv 
and 
eral. a 2A (> wa i 
V7A (r,t) 2s i ofA(r,t) = + We J(r,t) 
c? dt? 


which can be written in the following abbreviated manner. 


A(r,t) “Wy J(r,t) 

a Oe = =| 7 = 
(ITI-6) (v2 ~ = ) fee LE: 1... ye) 

ce? 3t? a = 
or 

a A(r,t) cy J (x,t) 
(ITI-7) Lior V(r,t) = J. (r,t) 
fy 


The solution to equation (III-7) can be written formally as 


follows: 
(III-8) eae = fa°r'* sat" Gir tice) “Wy cea 
V(r,t) : -1 aire") 
= 


where G(r,t|r}t') is the Green's function associated with the 
operator L, | and which must satisfy the following equation 


“a 


(III-9) Bet G(r,t|rjt') = 6(r-r') 6(t-t') 


where 6 is the Dirac delta function. 

In order to obtain a closed form solution to (III-8) one 
must find the form of the Green's function which is determined 
by (III-9). This can be easily done (see Appendix B) astaa 
Laplace Transform techniques to obtain the following solution 


for Gir tiric*} 3 


[eee ee a 
(III~10) G(z,t[rit') = - Zo St {t eo) 
|F - x | 
8 : ab 
The relation Hey. > oe has been used. 
c 
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Av othee requirement for a closed form solution to (III-8) is 

a functional form for J(r,t) and P (rit). The current source 
which is considered in this study is a single electron trav- 
eling with a velocity v(x,t) located at some point x at time t. 
Thus, Since this source iis localized rather than distributed 


through space the current density, as well as the charge 


density, can be represented using the Dirac delta function 


such that 

(III-11) g(rit') = ecB(r') &{r'-x(t") } ; 
and 

(ITI-12) py (ert) = ed{et-x(t')} 


where 8 = =, e is the amount of charge on an electron, and 
x(t) is the location of the electron at time t. Figure 8 
illustrates the geometry of the situation. Using the above 


expressions for J and p., and the solution of G(r,t[rjt") 


given in (III-10), the equations (III~8) for A and V become 


— een, ete FEET) 
(TII-13)  A(e,t)=-ye-/a°r sat! Se) B(x") é6{r'-x(t") } 
r-r’ 


and 


- ste HEE) 
(III-14) V(r,t)= er fatr'fdt' + a, 6{x'-x(t') } 
Vv r-xr' 


The equations given above could be solved #6 obtain a closed 
form for A and V. Then the operations indicated in equations 
(III-~2) could be performed to find expressions for E and B. 
However, it is less tedious if one simply performs the oper- 
ations mentioned directly on the integral forms of A and V 


given in (III-13) and (IIT-14). Performing these operations, 
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Figure 8. Geometry for a Single Electron Source. One should 
_ note that J(r,t) and On (r,t) are zero unless r(t) =x (t) 


one obtains, after much manipulation (see Appendix B), the 


following relations for E(r,t) and B(r,t): 


(III-15) E(r;t)= {| (n-B) (1- 8") 


ey 2 R2 ret 


nn. 


re | n a re 
+ a | OR x te) x Buleeyt | i 


anda 


(III-16) B(r,t) Lin x E(r,t)} 


where R=|R|=|r-r'| is the distance between the source (electron) 


a 


and the observer, n=5 is a unit vector indicating the direction 


between the source and the observer, 8 is the time derivative 


of the "velocity" of the electron (acceleration), and K=l-n-8. 
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The bracket with the subscript "ret" relates to the fact that 
the quantity in the brackets is to be evaluated at the re- 
tarded time, Pel iR EN 7el. That is, what the observer 
"sees" at time t is that which was created at an earlier time 
oa where the difference is the time that is required for the 
Signal to traverse the distance R. One should note that the 
fields indicated by (III-15) and (III-16) are divided into 
"velocity fields" which are independent of acceleration and 
"acceleration fields" which are acpenden’ on z in a linear 
fashion. One can show that the velocity fields are essen- 
tially sesere fields falling off as rR, whereas the 
acceleration fields are typically radiation fields falling 


off as RS 


(Ref 10, p 467). 

For long range effects of the EMP one is concerned 
with the radiated fields rather than the static fields 
Since the magnitude of the radiated field represents a flow 
of energy. Therefore, the fields of interest for the high 
altitude EMP depend on the acceleration of the Compton 
electrons as well as their velocity. That is, in order to 
solve for the radiated fields one needs analytical expressions 


for B{t). 


Recall that the primary source of the high frequency, 
high altitude EMP is the Compton electrons created by the 
gamma radiation from the weapon. Since the radiated fields 
depend on the acceleration of the electron, we shall consider 


two kinds of forces which can act upon the electrons 
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which are created in the absorption layer (30-50 Km). 

The first of these is the force which results from the inter- 
action of the electron (a charged particle) with the earth's 
magnetic field. The result of this interaction is a com- 
ponent of acceleration which is perpendicular to the 
instantaneous velocity of the electron. Secondly, the 
electron as it travels through the atmosphere, undergoes 

Many collisions with the surrounding atoms causing it to lose 
energy and consequently to de-accelerated. The result of 
these collisions is a component of neceleration which oa 
parallel? to the instantaneous velocity. Both of these 
forces will be explained in greater detail in the following 
paragraphs. 

Perpendicular Acceleration. The basic equation of 
motion for a charged particle under the influence of an 
electric or magnetic field is-the Lorentz force equation or 
(III-17) p = q(E + v x B) 
where p is the time dacivatige of the particle's momentum, 

q is the particle's charge, v is its velocity, and F and B 
represent the electric and magnetic field respectively. For 
the case of the electron in a magnetic field only, (III-17) 
_becomes 

(III-18) _ p= - e(v x B) (see Figure 9A) 
Using relativistic mechanics? ® one can easily prove the 


following identity between the momentum and the total 


9 ‘ 
For a zero order approximation 


lone should be aware of the fact that the electrons 
which produce the EMP are energetic enough to be considered 
relativistic 
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INCIDENT ELECTRON 


Figure 9. Electron Moving in a Magnetic Field. The (B) part 
of the above is the same as the (A) part except a 
definite set of coordinate axis has been defined. 


energy (€,) of a particle in a uniform magnetic field (Ref 15, 


p 28): 
| ee = & av 
p= &3 » at 
where &€; is assumed constant 
Using this relation in (III-18) one finds the following: 


(III-19) Ei dv 
, ae at 


If one defines the z-axis of a coordinate systen in the 


=- e(v x B) 


direction of the B field then (III~19) can be written in 


terms of the components of the vectors as: 


Be oe a er 
x f. oY 
. 2 
. ec*B 
(III-20) ser Ae aee rime 
Vo=0——> Ve _ (constant) 
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In order to solve the two coupled equations in (III-20) one 


2 
can let vee and then multiply the first equation by i=/-1 
1 
and add it to the second (Ref 15, p 59). This yields 
cms (Vo AV DS = 2, a) 
dt y x y x 


The solution to this is of the form: 


(III-21) V+ iv. =ge™ 
Y = - ~ia 
where g is a complex constant. One can write g as g=Ve = 


where ve and a, are real constants. Then using this in . 
(III-20) and separating real and imaginary parts one finds 
that 

(III-22) Vo = safe Sin(wt), MiG = V5 cos (wt) 

where a which represents the initial phase (t=0) has been 
taken as 0 for simplicity. From (III-21) it is obvious that 


the following relation holds: 


= sy cae — . . 
(III-23) Ne Wi + se |v| sin t (see Figure 9) 


Integrating (III-21) once more (Vv =9%) one finds that 


(ITI-24) x(t) = 5 aa r,cos (wt), y(t) = Yeo r, sin (wt) 
where 

Vv V8) 
(III-25) r= ~=- 

y ec’B 


Also if one integrates the equation for Vo in (III-20) the 
result if 


(ITI-26) z(t) a? + Vag ® 


Thus an electron moving in a uniform magnetic field with no 
other external forces will move in a helix (see Figure 10) 
having its axis along the direction of the magnetic field and 
with a radius ry given by (III-25). In the special case of 
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Figure 10. Trajectory of an Electron in a Uniform Magnetic 
Field. Part (B) is an “end on", cross-sectional 


. view of (A). 


the electron having no initial component of velocity in the 


direction of the B field (Vo,70), it moves along a circle in 


the plane perpendicular to the field (Ref 15, p 60). ‘The 


acceleration for the electron under these circumstances is 


also constant and can be found from (III~20) using (III-23) 


and (III-25). That is 


hak 


« 


or 


el 


— : vit = . a 
(11r-27)  [vj? = 2. = UV sin vw)" 


2 ; 2 


ere oy ae ae re 2) = 27 2 
Ve + Wy + Ve w (Vy + V_*) wv 


y oO 


: V & 


c*r,? 


1 se 


—* |B? = 


One can see that for an electron of a given energy the 


acceleration will be maximum when the incident electron is 


perpendicular to the B field (t=90°).. 


Thus the sin ~t term 


is seen to represent a scaling factor, the magnitude of 
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which depends on the angle of incidence of the electron 
with the magnetic field. 

Parallel Acceleration. When speaking of the "slowing 
down or "stopping" of an electron in a material such as air 
it is useful to define two quantities known as the stopping 
power and the range. The stopping power is defined as the 
amount of energy lost by an electron per unit length of 
path through the stopping material (Ref 6, p 185). The 
stopping power is symbolized by -d&/dx since it signifies 
the rate of energy change with respect to distance eevevea: 
The range of- the electron, or the length of the path it 


travels before coming to rest, can be defined in terms of 


the stopping power as follows: 


re) re) dt 
(IITI-28) Ry = f ax = - f at /ax 
ba Ey 
where R, is the range of an electron with initial energy of 


1 
—,. If one considers the case where an electron loses an 


increment of energy AE So: that its energy becomes &)=&,-AE, 
then one can compute the distance (Ry) that the electron 
must travel to lose this energy increment with the following 
relation: 
E2 
“1,2 7 7 a aE ax 


or separating the integration into two intervals 


fe) 18) 
Rp = Le dé }-(-Ss. a 


E1 d&/dx E2 d&/dx 


or using (III-28) 


ie 
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Rigo By ey 


(ITI-29) 

Equation (III-29) is necessary in this development 
because there are no simple, analytical expressions for the 
forces involved in the slowing down of an electron by a 
material. In fact, the majority of reliable information con- 
cerning this subject has been derived from experiment and is 
easily obtainable in the form of the total range (i.e. Ry or 
Ro and not Ri 2) of electrons with known initial energy. 
Therefore, in order to obtain information concerning the. 
acceleration of an electron in such a situation it will be 
necessary to construct a numerical model based on physical 
assumptions. 

The first assumption which shall be made in this model 
is that the electron remains on its original path while it 
is being slowed down and is not scattered in other directions. 
That is, it is assumed that the forces acting on the electron: 
are viscous in nature, changing only the energy of the elec- 
tron and not affecting its direction. This approximation is 
best for times soon after the creation of the Compton elec- 
tron and loses its validity after the electron has lost a 
large portion of its energy. As was mentioned, the forces 
‘acting upon the electron are assumed to be viscous in nature 
also implying that the slowing down procedure is continuous 
and that the acceleration function is a smooth curve. An 
empirical relation for the range of an electron as a function 


of its energy has been found using these assumptions. This 


relation is (Ref 17, p 16-4):. 
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2 
Joie. 


(III-30) Ree 
E, 40.26 


with + in Mev and Ry in gm/om? t+ 


In order to numerically approximate the eventos 
slowing down of the electron it will be assumed that the 
electron loses small increments of energy in a point-wise 
manner along its path. For example, if the electron has an 
initial energy of €,, and the energy loss increment is Ag, 
then it will be assumed that the electron travels a distance 
Ria! given by (III-29), with a constant energy f° After 
it has traveled this distance the electron then instanta- 
neously "loses" A& amount of energy thereby acquiring an 
energy of —,=€,-~Ac. This procedure is then repeated ina 
step-wise fashion throughout the slowing of the electron. 
Figure 11 gives a pictorial explanation of the situation. 
It should be pointed out that the PeGyenae of this approach 
increases as the magnitude of A& decreases. One should also 
note that the successive distances between energy increments 
(i.e. Ri 2° Ro 3! etc) are not necessarily equal. 

The acceleration can be approximated by giving the 


acceleration at the ith point (see Figure 11) such that 
(iets  ieit o (ae Aad 


(TII-31) Vv, = 
- ay = t., -¢ 
i+l ph | Led a=) 


x 


where V5 is the velocity corresponding to Ey which the 


linen measured in these units Ry is simply the actual 
distance traveled in cm multiplied by the density of the 
stopping material. This is convenient for experimental 


purposes. 
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Figure 11. Approximation of Continuous Slowing Down of Electron. 
The ith point is not at the point of energy change in 
order to avoid the question of infinite acceleration. 


particle has at the gh point, x5 is the total distance from 


th point, and t; is the time required 


the point i=0 to the i 
to traverse the distance Xi° If one realizes that the second 
term on the right hand side of (III-31) approximates the ve- 


locity which the electron has at the zth 


point, then the 

expression for ve can be written as: | 
: . Vv. -~ V. 
(III-32) Vv, a 
; i+]. 221 


_ {v;} 

' Therefore, one can use (III-32) for an electron of initial - 
Suergye by incrementing the energy as indicated above and 

using (III~29) and (III-30) for the values of x5 while sub- 
stituting the values for Vv; corresponding to each energy &;- 
This will yield a series of points for V from which a continuous 
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function can be implied. 

In order to illustrate the method described above, the 
calculations were performed for 1 MeV electron at an altitude 
of 40 Km. Figure 12 illustrates the form of the acceleration 
by showing the relative velocity.(8=V/c) as a function of 
time. Although the actual numbers are different for aeiSue 
altitudes, the general shape of all the curves is the same 
for the different altitudes. One should note that the major 
portion of the electron's acceleration occurs in a relatively 
short time at the end of the electron's trajectory. 

Electron trajectory. The electron in the atmosphere under- 
goes both the parallel and perpendicular acceleration simulta- 
neously. Thus, the developments of the previous two sections 
neea to be combined in some manner in order to give the elec- 
tron's true trajectory. Since it is not possible to derive a 
’ self-consistent simple force equation which would simulta- 
neously consider both effects, it is necessary to again create 
a numerical model based on physical assumptions. 

It will be assumed that the effect of the air molecules 
(pstaiidi eeesiecseion) on the electron is the same as given 
in the previous section. That is, the energy of the electron 
erenees by a definite increment at various points along its 
path. In between these changes, the electron moves at a con- 
stant energy and consequently at a constant speed. However, 
since there is a unieoRn magnetic field (perpendicular accel- 
eration) acting upon the electron, the electron does not 
traverse a straight line in between energy changes. Rather, 


the path forms a segment of a helix, the radius of curvature 
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of which is given by equation (III-25) since the speed of 
the electron (and energy) is constant over this path. The 
arc length or actual distance the electron travels for each 
segment between energy changes is given by equations (III-29) 
and (III-30). As will be shown later, an important variable 
is the angle which the instantaneous velocity vector of the 
electron makes with the initial direction (line-of-sight) of 


point 


the electron. If one measures this angle at the i 
(as defined in the previous section) then it can be given in 


terms of previously defined quantities as follows (see 


Figure 13): 
(III-33) y= = 
a 
z reo 
i=l 


th point. To give 


where vy is the angle as measured at the i 
the reader a picture of the electron's trajectory using this 
model, calculations were made for the various quantities 
defined at 40, 30, and 20 Kits? These results are summarized 
in Table 2. One should note the differences in ranges for 

the various altitudes. It should also be pointed out that the 
electron tends to spiral inward and at higher altitudes will 
make several complete revolutions before coming to rest (see 
wy. column in Table 2). 


1 


1eohese calculations assumed that the density of the air 


~“Z/Z 
/2., 


varies exponentially (i.e. p=. ). 
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Table 1. Calculations for Electron Trajectory. 
(|B|=10 * weber/m?) 


Electron Total Path . Total Time Radius of 
Energy Traveled Elapsed Curvature Va 
(MeV) (meters) (shakes) * (meter) (radians) 


ABSORPTION ALTITUDE = 40 KM 


1.00 0.0 0.0 47.49 0.0 

0.95 73.94 26.19 45.71 bead 
0.90 147.88 52008 43.92 3.17 
0.85 221.83 78.95 42.13 , 4,86 
0.80 295.77 105.55 40.32 6.61 
0.75 367.25 131.42 38.50 . 8.39 
0.70 441.19 O84 39 36.67 10.31 
0.65 512.67 184.67 34.82 12.26 
0.60 581.68 210.29 see95 14.24 
0.55 653.17 237.13 31.05 16.41 
0.50 419672 262.45 29.14 18.55 


ABSORPTION ALTITUDE = 30 KM 


1.00 0.0 0.0 47.48 0.0 

0.95 17.72 6.28 45.70 ‘o¥ 
0.30 35.44 12.58 43.92 -76 
0.85 53.16 18.92 42.13 1.16 
0.80 70.88 25.29 40.32 1.58 
0.75 88.01 31.50 38.50 2.01 
0.70 105.73 37.96 36.67 2.47 
0.65 122.86 44.26 34.82 2.94 
0.60 7139.40 50.40 23495 3.41 
0.55 156.53 56.83 31.05 34.93 
0.50 172.48 62.90 29.14 4.45 

ABSORPTION ALTITUDE = 20 KM 

1.00 0.0 0.0 47.48 0.0 

0.95 4.25 1.50 45.70 09 
0.90 8.49 3.02 43.92 18 
0.85 12.74 4.53 42.13 28 
0.80 16.99 6.06 40.32 236 
0.75 21.09 7.55 38.50 -48 
0.70 25.34 9.10 36.67 +09 
0.65 29.44 10.61 34.82 -70 
0.60 33.41 12.08 23495 82 
0.55 37.51 13.62 31.05 94 


0.50 41.34 15.07 29.14 1.07 


*1 shake=10 ° second 
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Figure 13. Electron Trajectory in Atmosphere. Part A is 
for the first two points. Part B is for the 


ha point and the hes point. 


. Fields for Electron Moving in-.Specific Manner 


The general. field equations for an accelerated electron 
have been derived and are given by equation (IITI-15). Further- 
more, it — been asserted that the electron is acted upon by 
a force which is perpendicular to its inet ancencous velocity 
(B - field) anda force which is parallel to this velocity 
(interaction with air rieleiles):, Since some insight has been 
gained onceeiing the nature of the acceleration of the elec- 
tron when subjected to these forces, it is now possible to | 


arrive at relations for the fields which will be radiated by 


_the electron as it is accelerated in this manner. 


First, however, it is useful to define a quantity known 


as the Poynting vector S which represents the energy which 
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flows across a unit area per unit time and is associated 
with radiated fields. if e and B. represent the radiated 


portions of the fields then S is defined as 


(ITI-34) S = a (E. x B.) 
V 
or using (III-16) 
=. 1 sic cay aace = 1 * me a: at F ees 
5 ci {E, x (n xX E,)} = cli, Sone E,) E,(n E,)} 


One can easily show that e is perpendicular to n which 


implies that 


(III-35) ° S& 


Using (III-35) one can then find the energy which is flowing 


in the direction of the observer by the following 


6 sy eat fe. |2 
CH, a 


or, uSing the radiation term in (III-15) 


_ * ee Ce “a PN tn a 
(I1I~36) {Sen}, = res) [5 In x ((n-B) x BHI?) 
Vv 


RoR" ret 
Using (III-36) one could then compute the amount of energy per 
unit area radiated during a finite period of acceleration, t' = 
' — 
Ty) and t T 
t=T 


2" and detected in a corresponding interval, 


+ {R(T,)/c} and t = T. + {R(T,)/c}, by the observer 


1 
with the following: 


2 


T,+R(T5)/c 7 
u= f AS niga ae 
T)+R(T))/c 
or | 
9 
(ILI-37) u=f (S+n) (st) at! 
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{ 
uineKe U is the energy measured. It is obvious from (III-37) 
that the quantity (en) (S,) represents the power (energy/time) 
radiated per unit area = the electron in terms of the "local 
time" of the electron rather than the retarded time defined 
earlier. Thus, with the dafinition of the solid angle in 
mind, one can define the power radiated per unit solid angle 
and measured by the observer at a distance R from the source 
as follows: 


dP (t') 


aq = RB’ (Sen) ( 


dt' 
or using the relation for retarded time (see Appendix B) 


aP(t') 
dQ 


which, using the definition of kx and equation (III-36), 


= KR2 (Sen) 


becomes 


a Nn ° 2 
aP(t') _ 1 ee,” ;|n x {(n-8) x B}] 
(III-38) -—3.— = -— (+— { ; x a I 
an ce, ia iy a. B)° 


One can reduce (TII-38) to a simpler form for the case 
corresponding to an acceleration which is acting parallel to 


the instantaneous velocity by recognizing that 


e tad a e 


(III-39) ‘ x {(n-8) x B} = . #inxb-oee)= n x (nxp 


Since 8 is parallel to 8. Using (III-39) and letting w be 


the angle between 8, 8 and n (see Figure 14 A) then one 


finds that 
U ° 
(ry39)). Sr ce (28)? _sin? 
an - ec 4t (1-8 cos wv?) 


where Pe signifies the power radiated which is due only to 


parallel acceleration. A similar reduction in form is 


possible for the case where the acceleration is perpendicular 
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j 
to the instantaneous velocity. If one defines n as the angle 
between 8 and the projection of n into the plane perpendicular 


to 6 (see Figure 14 B) then one finds that 


(III-40) dP (t") Pee (28° 1 {1 sin? wv cos? yn 
fey ee eee ~ ee eee Ne: 
dQ ce, 4t (1-8 cos tb) y2(1+B cos )2 


where y?=(1-B2) 7" and Po Signifies the power radiated which is 
due only to perpendicular acceleration. 

With regard to the high altitude EMP the next step is to 
investigate the relative eentraipations of radiation due to 
parallel acceleration and perpendicular acceleration. Cal- 
culations can be made for (III-39) and (TII-40) at various 
points along the electron's trajectory using the model out- 
lined in the last part of the previous section to obtain 
values for the acceleration terms at various points. If one 
assumes that the observer is situated on the original line-of- 
sight, then n=180° and wp at each point is the angle vs defined 
in the previous section. One should note that n=0° or 180° 
corresponds to the minimum value for (ITI-40) with respect 
to n. The calculations indicated here were performed by the 
author for Compton electrons created at various altitudes 
ranging from 60 to 20 Km. The results at 40 Km for early 
times are presented in Table 2 on Page .49. As one can 
readily see, the radiation resulting because of the electron's 
interaction with the magnetic field (AP ./d&) far outweighs 
that resulting from the linear slowing down of the electron 
(AP, /d2) . Although Table 2 shows only the calculations at 


early times, and at one altitude, the results are found 
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Figure 14. Geometry for Angular Distribution of Electron _ 
Radiation. Part A is for parallel acceleration; 
Part B concerns perpendicular acceleration. 

to be the same at late times and at all altitudes. Thus, 

based on these calculations, one Gan. Gag that the electro- 

magnetic radiation which constitutes the high mireteaae EMP 

results from the turning of haghh energy Compton electrons by 

the earth's magnetic field. oo 

It is beneficial to investigate even further the nature 

of the radiation resulting from the turning of the electron. 

Table 2 snows that for a period over which the radiated power 

drops four orders of magnitude, the velocity, 8, of the elec- 

tron is essentially constant to two significant figures. Thus, 

although much time has been devoted to the slowing of the elec- 

tron in (Hel-atmesphare’”; it is seen here that, with regard 


l3nnis development will be very important later, however, 
when discussing air conductivity. 
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Table 2 Comparison of 


Electron 
Energy 


(MeV) 
1.0000 
-9964 
~9929 
- 9893 
~9857 
-9821 
~9786 
- .9750 
.9714 
-9679 
~9643 
- 9607 
~9571 
- 9536 
“.9500 
~9464 
~9429 
9393 
-9357 
9321 
~9286 
-9250 
-9214 
-9179 
.9143 
~9071 


B 
(=) 
~ 94) 
-941 
- 940 
-940 


940 


940 
£939 
.939 
.939 
.938 
938 
938 
937 
937 
.937 
936 
936 
936 
935 
.935 
.935 
934 
.934 
934 
.933 

2933 


dP 


({Bl=1x10°° webers/m?) 


Total Time 


Elapsed 
(shakes) 

0.00 

1.82 

3.64 

5.46 

7528 

9.10 
10.92 
12.75 
14.57 
16.39 
18.21 
20.03 
21.85 
23.67 
25.49 
27.32 
29.14 
30.96 
32.78 
34.60 
36.43 
38.25 
40.07 
41.89 
45.54 
47.36 


dan 


dP 


and sa for Electron at 40 Km. 


AL 


wv. 


a 
(radians) 
0.000 
-108 
sad 
2325 
434 
~544 
+053 
- 763 
-873 
984 
1.094 
1.205 
1.316 
1.428 
1.539 
1.651 
1.764 
1.876 
1.989 
2.102 
2.216 
24329 
2.443 
2.558 
2.787 
2.902 


(4) 
0 
1. 758-39 
2.24E-39 
1.17E-39 
4.56E-40 
1.69E-40 
6.46E-41 
2.62E-41 
1.14E-41 
5.27E-42 
2.59E-42 
1.34E-42 
7.30E-43 
4.14E-43 
2.43E-43 
1.47E-43 
9.10E-44 
5.77E-44 
3.73E-44 
2.41E~44 
1.58B-44 
1.03E-44 
6.58E-45 
4.10E-45 
1.29E-45 
5.21E-45 


ap, 
da 
(*) 
1.07E-13 
9.32E-14 
4.81E-14 
7.00E-15 
6.80E-17 
2.98E-16 
4.74E-16 
3.89R-16 
2.69E-16 
1.762-16 
1158-16 
7.64E-17 
3.638217 
2.61E-17 
1.93E-17 
1.475-17 
1.148-17 
1.10E-18 
7.42E-18 
6.20E-18 
5. 28E-18 
4.60E~18 
4.08E-18 
3.70E-18 
3.21E-18 
3.08E-18 


*These calculations were performed using CGS-Gaussian units 
rather than MKS units. 
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to the high altitude EMP, the radiating electron can be con- 
sidered as moving with a constant speed and being turned by 
the magnetic field of the earth for the time periods of inter- 
est. Making this assumption, one can then calculate the angle 
that the instantaneous velocity vector makes with the initial 
direction of the electron (the "line-of-sight" direction for 


example) at any time t as follows: 


(111-41) w(t) = BSE 
xr 
B 
where ) is the indicated angle and r is the radius of curva- 


B 


ture given by (III-25) corresponding to an electron with 

speed 8. Using (III-41) in (III-40) one could obtain in- 
formation concerning the power or the magnitude of the E 
fields which an observer located along the initial direction 
of the electron would measure as a function of time. Such 

a calculation is illustrated in Figure 15 which is a graph 

of normalized values of the time dependence indicated by 
(III-40) and (III-41). An even more intuitive understanding 
of the distribution of this radiation is gained from Figure 16 
which illustrates an iso-field line for an electron with a 
certain instantaneous velocity. From these two figures, one 
can see that a large portion of the electron's radiated en- 
ergy is directed forward in the direction of its instantaneous 
velocity. Also, it is evident that the strength of the fields 
drops off rapidly as one moves away from this direction. thus, 
an observer measures an initial maximum field strength which 
drops off drastically after the electron turns only a few 


degrees. 
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FIG.15, E-FIELD FOR SINGLE ELECTRON 
MEASURED BY LINE“OF-SIGHT OBSERVER 
(NORMAL IZED). 
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Line-of-Sight Effect. Since there are many Compton 
aieetrons created in the upper atmosphere by a high altitude 
nuclear detonation, it is necessary to know which of these 
contribute to the high frequency EMP measured by an observer 
at a particular location. It has been previously asserted that 
the Compton electrons which contribute to the high frequency 
(early time) EMP are those which are created by gamma photons 
moving on the "line-of-sight" direction for the observer.with 
respect to the weapon detonation (see Figure Vs tas There are 
two primary reasons for this which involve the time phasing 
of the radiated fields and the angular distributions of the 
radiated energy. . 

The important factor with regard to the time phasing 
aspect is that the Compton electrons are traveling at essen- 
tially the speed of light aaving the time they radiate the 
Major portion of their energy. The importance of this fact 
can be understood by considering two gamma photons emitted 
at the same time and moving together on the line-of-sight, 
each of which creates a Compton electron but at different 
times (see Figure 17). Since the first electron created is 
essentially moving at the same speed as the remaining gamma 
photon, the fields which are radiated in the time interval 
T (T, small enough so that the electron has not appreciably 


1 
turned off the line-of-sight) and which propagate along the: 


Mone reader is reminded that it has been assumed that 
all Compton electrons initially move in the incident direction 
of the gamma photon. 
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Figure 16..: Radiation Pattern for Electton in Circular Motion. 
(Note that the coordinate system which defines v 
(and n) is constantly changing since the z- axis 
is defined by B(t).) 
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line-of-sight constitute a pulse the "leading edge" of which 
is almost even with the "leading edge" of the gamma photon. 
Similarily, Oe ere the electron created later by the remaining 
gamma photon is also moving with approximately the speed of 
light, the fields which it radiates in the time interval Ty 
and which propogate along the line-of-sight constitute a 

pulse the "leading edge" of which is almost perfectly matched 
up to the first pulse. These two pulses, then, will arrive at 
the observer in such a manner as to result in almost perfect 
constructive interference. One can infer from this argument 
that the more the electron's instantaneous velocity vector 
deviates from the line-of-sight the less the component of 
velocity in the line-of~sight direction becomes and the longer 
the distance which the radiated energy must travel to the ob- 
server becomes. The result of this is that the contribution 
of the radiated pulse to that which the observer measures 
becomes less and less significant as compared to the con- 
tribution of the line-of-sight siecevens until it is nothing 
more than background. Thus, when speaking of "line-of-sight" 
electrons with respect to time phasing one included electrons 
created within a small cone centered about the line-of-sight 
such that their valocity component parallel to the line-of- 
sight is approximately equal to their total velocity and 

such that the distance which the radiated energy travels to 
the observer is approximately equal to that distance which 
the energy radiated by actual line-of-sight electrons travels. 


Independent of the time phasing and more intuitively 
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obvious is the effect of the angular distribution of the 
radiated fields. From Figuve 15 and Figure 16 one can see 
that the electron which contributes the most to the power 

the observer measures is the one which is moving directly 
toward the observer. Referring to Table 2 one also can see 
that when the angle between the instantaneous velocity and 
the direction to the observer is just a few degrees, the 
radiated fields measured by the observer have dropped by 

an order of magnitude. Thus, the electron which is created 
by a gamma whose incident direction is just a few degrees 
from the line-of-sight direction will not significantly 
contribute to the early time EMP. It can be argued that an 
electron with initial direction away from the observer can 

be turned by the earth's magnetic field such that is eventually 
"points" at the observer. This argument is valid insofar as 
the motion of the electron is such that it satisfies the time 
phasing arguments given in the previous paragraph. 

Although it has been shown that the important contribution 
to the high altitude, high frequency EMP is made by line-of- 
sight Scéeeoas: one should realize that, when speaking of the 
electrons which contribute to the fields measured by an ob- 
server at a particular location, some electrons which are 
created by gamma photons with directions slightly off the 
line-of-sight are to be included. The degree of deviation 
from the line-of-sight is determined by the constraints men- 
tioned in the two previous sections. It is these qualifications 


which should be kept in mind when speaking of the high altitude, 
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high Eeaquency EMP as a "line-of-sight" phenomenon. 

Super imposing Effect. When computing the net effect of 
all the electrons weet contribute to the EMP one is primarily 
concerned with the total field strength which the observer 
measures as a function of time. The variation in time of the 
field strength will depend, among other things, on the time 
nature of the gamma photon output. From equation (II-10) in 
Chapter II it has been shown that the number of gamma photons 
passing through a unit area at some distance R from the burst 


and at some time t" is given by 


E> 
(ILI-42) ®(R,t") = _ £ (T") F(R) 
Y 
where ~HOB/z . R/z, 
ie {Ze (l-e )} 
41R? 
and 
ne R 
c 


Furthermore, using (II-16) one can write the number of Compton 
electrons which are created per unit volume per unit time at 
time t" as 


dN. (8, t") Y 
pri = ory. " 
ab EY A(R) F(R) £(T") 


The number of these electrons which would be created in a 
time interval dt can be found from the following: 
Y 
(III-43) an. (R,t") = =< A(R) F(R) E(T")at 
pri EB 
y 
Tf EB, (rst) represents the general time dependence of the 


field strengths radiated by a single electron and measured by 
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an observer at a distance r away and at some time ae then 


for an electron created at t" the fields measured at time t 
after the weapon detonation would be given by E,(r,t-t"). 
Thus an observer would measure the following for (VedN ys) 
electrons created a distance r away from him and R from the 


burst: 


E, (r,t) = VeGN 3 (ReETVE, (rs t-t") 


where V is the volume enclosing the electrons which contribute 
as earlier described. If the time origin is shifted such that 
t=0 when T"=0, then the above can be written as 


B, (ret) zs VAN: EMRE JE (x, t-T ) 


or using (III-43) and recognizing that the r dependence for 
E, is l/r for radiated fields, 
VY 


: eis ile Mee Aeteaea 
(EII-44) | E, (r,t) rE, A(R) F(R) {E(T")E, (t-T") aT" } 


E, (r,t) just accounts for those electrons created within an 
infinitesimal time increment of dT. There will be electrons 
created before and after this increment which will contribute 
to the total field strength, E, (x,t), measured by the ob- 
server. Thus, in order to obtain the total field strength 
caused by Compton electrons created within a finite volume! ® 


centered about a point some distance R from the burst and 


Measured at some time t by an observer located a distance r 


lSohese fields were created at t'<t_-= (retarded time) 


16por further information concerning this volume see 
Chapter IV. 
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from the electrons, one should integrate (III-44) from 


T=0 to T=t-2. That is, 


r 
VY_A(R)F(R) ‘"e 
(III-45)  E, (r,t) So f £(@")B (t-T")arT" 
y~ oO e 


As has been indicated, equation (III-45) is for the 
electrons created within a certain volume centered about a 
point some distance R from the burst. One can sum over all 
such points by letting V=A*dr where A is the epsaaraeceioust 
area of V lying ina plane perpendicular to the line-of-, 


a and dR is an infinitesimal length parallel to the 


sight 
line-of-sight, then integrating (III-45) over the absorption 
layer. Letting A=A(R) one can then write the following for 


the fields which an observer located a distance X from the 


burst would measure at time t. 


(X-R) 
: YX A(R)A(R) F(R) oe 
(III-46) “E(x,t) = gt s aR ay fate (2) E, (t-7") 
yo Oo 


One should note that it was not necessary to make a correction 
for time lags between altitude because of the time phasing 


which was mentioned previously. 


lone constraints mentioned in the previous section 
regarding the line-of-sight effect determine the size of A. 
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Iv. Sample Calculations 


The purpose of this. chapter is to use the theory devel- 
oped in the previous chapters and obtain data regarding the 
nature of the high altitude, high frequency EMP. The tech- 
niques used to solve the various equations and arrive at the 
physical quantities mentioned previously will also be demon- 
strated. For simplicity, all these calculations involve an 
observer directly beneath the burst and who is located at 
the magnetic equator. Thus the incident velocity of the- 
Ganseon electrons will be perpendicular to the earth's 
magnetic field. Furthermore, it shall be assumed that pie 
earth's magnetic field is uniform over the absorption layer 


and has a field strength of 1.0x10 * webers/meters. 


Line-of-Sight Area 

As was mentioned in the previous chapter, the high 
frequency EMP is “almost" a line-of~sight phenomenon. That 
is, there are electrons just slightly off the line-of-sight 
which contribute significantly to the EMP measured by an 
observer. The first task in computing the EMP, then is to 
determine the number of the electrons which contribute to 
the observer as a function of distance from the burst. For 
these calculations, it shall be assumed that all the elec- 
trons which contribute, do so in exactly the same manner 
as the electron which is precisely on the line-of-sight. 
Thus, the cross-sectional area located at a certain distance 


R from the burst which contains all the electrons that con- 


59 


GNE/PH/74-7 
i 


bh epdee will be considered to act as a point source to the 
observer. With regard to this approximation, it is obvious 
that the effect of electrons which originate with a velocity 
slightly off the line-of-sight yet which are turned so as to 
"point" at the observer wit be underestimated. Similarly, 
the effect of the electrons originating with a velocity off 
the line-of-sight and which are turned away from the observer 
will be overestimated (see Figure 18). To obtain a measure 
of the error introduced by this, equation (III-40) with n=180° 
was integrated bwtween an angle be and Yo + 6° for ve ranging 
between -1° and 1°. The per cent difference between each 
integration and the integration between We 0%. and Yo + 6° 

was then computed. These results are presented in Table 3. 
For the purpose of this study, it was decided that the error 
introduced by be = + 0.5 was tolerable. Thus, the maximum 
angle between the incident direction of the electron and the . 
line to the observer in order for the electron to be con- 
sidered as a line-of-sight electron will be set at +0.5°. 

If one assumes that the cross-section area A(R) at each 
altitude is a circle, then one must find the radius of this 
circle in terms of W, (=40.5°) at each altitude. Using the 
. geometry indicated in Figure 19, with the observer at sea 


level and letting X be the sought after radius, one finds 


LOS 
that 


. _ _ HOB / > HOB \2 
(Iv-1) Xiog (R) = 2tany, ” Stan.) PRE 
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Table 3 
Difference of Total Power 


Received for Electrons 


off Line-of-Sight 


Ye % Yo 
| DEGREES DIFF DEGREES 

~1.0 1132 O.1. 
-0.9 10.3 0.2 
-0.8 9.3 0.3 
-0.7 8.2 0.4 
-0.6 de? 0.5 
-0.5 6.1 0.6 
-0.4 4.9 0.7 
-0.3 3.8 0.8 
-0.2 2.5 0.9 
-0.1 1.3 Lae 
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(4) 0909 
B 

ELECTRONS WILL 

TURN TO THE RIGHT 


IN This as 
te Yn6 0 AND Fes 


Pi] OBSERVER 


} 


Figure 18. Iso-Power Lines for Electrons in Three Different 
Configurations Relative to Observer. Ina Set. 
time period during which each electron turns the 
“same distance and consequently sweeps out the same 
number of degrees (at least >Prg?) then the con- 


tribution of A will be greater than B which will 
be greater than C. 


Therefore, for a or ecular ences sectional area A(R) one has, 


: 2 
nin) = te (HOR SREY me 
; Oo 


A graph was made for xX as a function of altitude for a. 


LOS 
weapon detonation occurring at 100 Km and with Wo=0.5°. This 
is given is Figure 20. : 

The assumption that the electrons enclosea in the circular 
area A(R) above the observer and centered. on the line-of- tea 
act as a point source lying on the line-of-sight ignores the 


n dependence in equation (III-40). For the electron orig- 


inating on the line-of-sight, n is always equal to 0° or 180° 
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HOB=ER+Z 


OO 0 
B - 
(ouT OF PAPER 


Figure 19. Geomerty for Finding the Line-of-Sight Radius. 
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So FIG.20, LENE-OF-SIGHT RADIUS AS 

o A FUNCTION OF ALTITUDE, 
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and cos*n=l. However, there are electrons included in A(R) 
which have values of n ranging from 0° to +180° (see Figure 
21). By saying that the disk with area A(R) appears as a 
point to the observer, one is implying that the distance 
Xhos cannot be revolved. Therefore, as far as the observer 
is concerned, the electrons all appéar at a point and have 

a value of n=0° or 180° depending on the orientation of the 
magnetic field. To test the validity of this assumption one 


can compare the magnitude of X og (8) at various altitudes to 


L . 

the magnitude of the distance to the observer. For example, 
using Figure 20 one can see that at an altitude of 50 Km the 
value of Xios is about 0.22 Km. The ratio between this and 
the distance to the observer is .004 so that one could say 
that Xnos<<2 where z, in this case, is the distance to the 
observer. Similarly, at 20 Km the ratio is .007. Therefore, 
the assumption of a point source is a relatively good one. 
One should note that the A(R) which meets the point source 
specifications falls well within the line-of-sight require- 


ments mentioned in the previous chapter. 


Integration Over Absorption Layer 

The integration over R in equation (III-46) can be 
approximated by summing over segments of the absorption layer 
each centered about some distance Ry and having a depth of 
AR. The quantity AR must be such that the macroscopic cross 
sections and all other quantities depending on R are constant 
over the interval, since all quantities in this interval will 
be assumed to have the value given at R,- It was found that 
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EE eee eC 
Figure 21. Top View of Electron in Line-of-Sight Area A(R). 


AR=1 Km was a relatively good approximation. Using this ane 
can now compute the magnitude of the fields resulting from a 
‘single electron at R, which an Siacyer would measure at a 
certain time, multiply that number by the Compton electrons 
per unit volume formed at R, at an epeeuiake earlier time, 
anda then multiply by A(R,) AR to get the total effect of (a 
-eylindrical volume of electrons centered at R,. This pro- 
cedure could then be pepearsd for all other Ry in the absorp- 
tion aves and the various contributions summed up in order 


to obtain the net effect. 


Discussion of Results 
The calculations which have been previously discussed 
were performed for a nuclear detonation occurring at.100 Km 


with a gamma yiéld of one kiloton. The magnitude of the E 
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fields which an observer, located directly beneath the burst 
at sea level and at the magnetic equator, would measure, based 
on these calculations, is presented in Figure 22. The length 
of the pulse is about 5 shakes (1 shake = 10 * sec) and reaches 
peak fields of the order of 12.0 x 10"* volts/meter. The 
validity of the section of the graph for times after ten 
shakes is doubtful since by this time many of the electrons 
have slowed considerably below the value of 8=0.94 which was 
assumed for the calculation. Ina paper by Karzas and Latter 
(Ref 11) using a macroscopic approach with similar weapon 
parameters, peak field values of 6.0 x 10" volts/meter are 
reported. Thus, the microscopic approach does give answers 
which are in the same order of magnitude as accepted theories. 
The difference between the two could be accounted for by the 
fact that conductivity has been ignored here but was accounted 
for in the Karzas and Latter paper. 

Although the previous calculations were done for the case 
where the electron's incident velocity is perpendicular to the 
earth's magnetic field, it is a simple matter to show the ef- 
fects of different orientations using this microscopic model. 
Equations (III-35) and (III-40) indicate that the magnitude of 
the E fields measured by an observer are proportional to the 
acceleration, 8. From equation (III-27) one can see that for 
an electron of a given energy the acceleration is proportional 
to sin*t where t is the angle between the incident direction 
and the magnetic field. Thus, the magnitude of the fields 


measured by an observer will vary as sin*t according to his 
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FIG,22,. EMP MEASURED AT GROUND AS 
A FUNCTION OF TIME, 
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eaaseie bearing with respect to the weapon burst. This 
point is illustrated by Figure 23 on the next page. The ob- 
server at point E will measure the maximum field magnitudes 
since the angle, Tt, between the initial direction of the 
electrons and the magnetic field is 90°. However, the ob- 
server at point A will not measure any fields since the 
direction of the electrons' motion is parallel to the magnetic 
field (t=0°). - The field strengths measured by the observers 
at points B, C and D increase as sin*t according to the 
respective values of t in each case. 

One should be aware of the potential range of the high 
altitude EMP which can be inferred from its "line-of-sight" 
nature and the facts of the previous paragraph. Since the 
high altitude, high frequency EMP is essentially a line-of- 
sight phenomenon, the effects of a burst located at a proper 
location with respect to the earth's magnetic field could be 
dstested from "horizon to horizon". This fact is illustrated 
in Figure 24. As one can see from this figure, the size of 
the Paes in which the EMP could be detected, for a given 
orientation with respect to the earth's magnetic field, 
depends on the height of the weapon burst. 

Calculations similar to that mentioned at the beginning 
of this section were performed for various gamma yields and 
heights of burst for the weapon. These calculations also 
ignored conductivity and considered only the gamma induced 
(ignoring x-ray induced) EMP. The results of these are 


presented on the following pages in graph form. 
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WEAPON BURST - 


| OBSERVERS ; | 


Figure 23. Effect of the Earth's Magnetic Field Orientation 
on the High Altitude EMP. 


Figure 24. The Horizon-to-Horizon Nature of the EMP. 
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V. Effects of Air Conductivity 


In the previous chapters it was assumed that the con- 
ductivity of the air surrounding the radiating electron could 
be neglected. The purpose of this chapter is to study the 
validity of this assumption. This will be done by briefly 
studying the physics involved in air conductivity with re- 
spect to the mechanisms relating to the creation of the 
conductivity and its subsequent buildup in time. There will 
also be a short section included to indicate the equations 
which must be solved in order to obtain the fields radiated 


within a conducting media. 


Formation and Buildup of Conductivity 
The basic cause of conductivity in the air is the 

ionization of the molecules of the air caused by the high 
energy Compton electrons. The result of this ionization is 
to release large numbers of low energy electrons (secondary 
electrons) into the media. Although these electrons are not 
energetic enough to add appreciably to the radiated fields, 
they do cause the environment to act as a conductor since 
they are free electrons. One can express the conductivity 


of the air in terms of the number density of these secondary 


electrons as follows (Ref 13, p 101): 


(V-1) o(R,t) 


where o is the conductivity measured at a distance R from the 


weapon detonation at time t, N, is the secondary electron 


ec 
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Hone e is the magnitude of the charge on a single electron, 
Mm, is the mass of an electron, and Vo is the frequency of 
collisions at an altitude z above the ground. It should be 
mentioned that, in general, the conductivity depends on the 
density of negative and agaitive ions as well as free elec- 
trons. However, in this case the contribution of the ions is 
negligible as compared to that of the electrons. The reason 
for this is that the conductivity is proportional to the sum 
of the inverses of the mass of the electron, positive ions, 
-and negative ions. Since the ions are orders of magnitude 
more massive than the electron, the terms involving the 
inverse of their mass will be negligible compared to the 

term involving the inverse of the electron mass. 

If one assumes that ve is constant with respect to time, 
then the only time dependence in (V-1) is that in the Nags 
term. There. are several mechanisms which affect the time 
dependency of the density of the secondary electrons. One 
such mechanism is, of course, the rate at which energy is 
deposited by the Compton electrons to ionize the air, there- 
by creating the secondary electrons. However, the chemistry 
of the air is such that there are also mechanisms whereby 
these secondary electrons may be removed. The mechanisms 
include the attachment of a free secondary electron to an 


0. atom to form an 0.- ion and the recombination of an 


2 


05+ ion and a Not ion with a free electron. Since the main 
purpose of this chapter is to test the validity of ignoring 
the conductivity, the removal mechanisms for secondary elec- 


trons will be ignored in order to arrive at a "worst possible 
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case" for conductivity. Thus, a lower limit for the time 
during which conductivity can be ignored will be developed. 

In order to pete ans number of secondary electrons 
created to a given amount of energy deposited by the Compton 
electrons, one can geting an average ionization energy W 
which is defined as the ratio of the energy lost by the 
Compton electron to the total ionization produced by it. That 
is, W represents the average energy required to be deposited 
in order to create a secondary electron. An accepted value 
for W in air is 34 eV (Ref 19, p 5-6). If a Compton electron 


was created at some time T and if & (ti) represents the 


DEP 
amount of energy which this electron has deposited at some 
time te after its creation, then the number of secondary 


electrons at time t is given by 


Epgp (t-T) 


(V-2) Nec (t) = W 
It is obvious from this argument that the number density of 
secondary electrons at a distance R from the burst at a given 
time t will be dependent upon the number density of the Compton 
electrons at this point in space and time. That is, at any 
time t after detonation the number of secondary electrons at 

a certain R will depend on the number of all the Compton 
electrons which have been created peevicus to t and the amount 
of energy which each has deposited between the time of its | 
creation and t. If ape represents the rate of change 


of the primary electron number density at R and T, then one 


can express the secondary electron number density at any time 


rae 
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t as follows: 


t aN. .(R,t) 


: gE (RET) 
_ = pri DEP 
(v=) NU (Rpt) = : { } }ar 


ee W 
where an R dependence has been included in DEP to account 
for any variation with altitude. 

In order to obtain numerical values for (V-1) one needs 
an expression for Vo (Ze The expression for Vag can be shown 
to vary with altitude primarily as the density of the air 
varies with altitude. Thus, if one assumes that the atmos- 
pheric density varies exponentially with altitude then one 


can write: 


(v-4) Vv, (4) = vee 2/20 


where vs is the collision frequency at sea level and Zo is 
the scale factor introduced in Chapter II. Using geometry 
indicated in Chapter II one can write the following for 


(v-4) 


vo (2) =iyuen HOBRRIZs 


Using (V-1), (V-3), and (V~-4) calculations were made for 
the conductivity as a function of time at various altitudes 
and for weapons with different amounts of gamma yields. A 


functional dependence for € (R,t) was obtained at each 


DEP 
altitude from data such as that indicated in Table 3. The | 
data was "fit" using a least squares technique and the re- 
sulting functional form was used in the integral in (V-3). 
The results of these calculations for an altitude of 40 Km 
and varying gamma yields are presented in Table 4 on the 
following page. To illustrate the variation in altitude, 


77 
% 


GNE/PH/74-7. 


Time 


(shakes 


e  e# 8 © * @ 


. . . 


CUMS AHAHUUBRWNNE © 
ON PHAWON HL HDWON & 


= 


« e e . e e ry . e 


OWDWAHONPWNYNE OO 


. e 8 «€ 


ONADAON SE AWOND 


bw 
° * e ° 
ON FPHDDON HL HADON S 


e ss oe @ 


+ * 


~ 
COW WDANAAHUAWNHDNE SO 


*y=2.0x10— 


) 


7 


Conductivity as a Function of Time 


sec 


~1 


Table 4. 


(Altitude 


SDEP 
(eV) 


GAMMA YIELD 


0 
1.403E02 
1.865E03 
3.586E03 
5.302E03 
7.014E03 
8.721E03 
1.042804 
1.212504 
1.382E04 
1.551E04 
1.720804 
1.888E04 


GAMMA YIELD 


0 
1.403E02 
1.865E03 
3.586E03 
5.302E03 
7.014E03 
8.721E03 
1.042E04 
1.212E04 
1.382E04 
1.551E04 
1.720E04 
1.888E04 


GAMMA YIELD 


0 
1.403E02 
1.865E03 
3.586503 
5.302E03 
7.014F03 
8.721E03 
1.042E04 
1.212E04 
1.382E04 
1.551E04 
1.720E04 
1.888E04 
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40 Km) 


ce) 


(Mho/Meter 


1 KILOTON 

0 

0 
5.314E-09 
1.082E-07 
2.303E-07 
3.627E-07 
5.010E-07 
6.427E-07 
7.867E-07 
9.320E-07 
1.078E-06 
1.225E-06 
1.372E-06 


5 KILOTONS 

0 

0 
2.657E-08 
5.410E-07 
1.152E-06 
1.814E-06 
2.505E-06 
3.214E-06 
3.934E-06 
4.660E-06 
5.390E-06 
6.125E-06 
6.860E-06 


10 KILOTONS 


0 

0 
5.314E-08 
1.082E-06 
2.303E-06 
3.627E-06 
5.010E-06 
6.427E-06 
7.867E-06 
9.320E-06 
1.078E-05 
1.225E-05 
1.372E-05 


3.002E-05 
6.115E-04 


1.301E-03 


2,049E-03 
2.830E-03 
3.631E-03 
4.445E-03 
5.265E-03 
6.091E-03 
6.920E-03 
7.750E-03 


0 

0 
1.501E-04 
3.058E-03 
6.505E-03 
1.025E-02 
1.415E-02 
1.816E-02 
2.223E-02 
2.633E~02 
3.046E-02 
3.460E-02 
3.875E-02 


0 


0 
3.002E-04 
6.115E-03 
1.301E-02 
2.049E-02 
2<G30B=-02 
3.631E-02 
4,445E-02 
5.265E-02 
6.091E-02 
6.920E-02 
7.750E-02 
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the conductivity at a given instant of time is plotted 


- +» vas a function of altitude in Figure 28. 


"Relative Significance of Conductivity 
To gain a deeper understanding for the significance of 


the values listed in Table 5, it is necessary to consider 


regard to the equation which involves the current density, 


(v-5) VxBeureh 
Vv 2 


one can express the current density as the sum of the Compton 
current density and the current density of the secondary 


electrons as follows: 


J=9 + J 
comp sec 


| In the previous chapters, it has been assumed that Tent: 
If one assumes that Ohm's law applies to the secondary elec- 


tron current density, namely that 


then (V-5) can be written as: 


1 2%E 


(V-6) VxBe=uJ + LOE + > BE 


v- comp 
Using (V-6) and the remainder of Maxwell's equations as 


indicated in Chapter. III and Appendix B, one can obtain 


potential V such that:78 
2 a < Wet eas 
b oa ee oe Ryo R/_| | 
Bede SRE eecgat ee rtd 
; e, c 
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Assuming that o is constant with respect to space and 
time. ; 
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If one applies a Fourier transform with respect to time to 
the left hand side of (V-7) the time derivatives can be done 
explicitly to yield the following for the operator L enclosed 
in the brackets: 

tan a oe ae _ iw)? 
(V-8) a ae iw yo } 


FT a? 


where Lam indicates that the Fourier transform has been applied, 
izV¥-1, and w» is the transformation variable. It can be shown 
that w corresponds physically to the frequency of the radiated 


= 2 
fields. If one divides (V-8) by fiw)" the results are: 
c 


OF ON: «lay ee is eo . 
Gio)? ‘rr (iw) 2 ve i) oa 


where eral/i eu has been used. Thus, one can obtain an esti-~ 
Mation of the relative importance of the term involving 


conductivity by comparing the magnitude of the term — 
Vv 


to 1. Since the actual fields which are radiated form a pulse 
which consists of many frequencies, it is necessary to arrive 
at some "characteristic frequency” of the pulse in order to 
obtain numbers for the aT term for various values of o. A 
erude approximation diien, is consistent with the Uncertainty 
Principle is to let w be equal to the inverse time width of 
the radiated pulse. In Chapter IV it was shown that the 


pulse was about 5x10 ° 


seconds long. Thus, w=2xl10’ inverse 
seconds. The results of such calculations are indicated in 
the last column of Table 4. If the number .01 is arbitrarily 


set as the level at which an is to be considered significant, 
Vv 


then one can see from Table 4 that the conductivity can be 
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ignored for a gamma yield of one kiloton during the times 
which the high frequency EMP is being radiated. However, as 
the gamma yield is increased to five and then ten kilotons, 
the conductivity does become significant at relatively early 
times. Therefore, it is fairly safe to ignore the conduc- 
tivity for weapons with relatively low gamma yields (< 1 KT). 
One should be aware that this last statement is based on 


ignoring .01 as compared to l. 


- 


Field Equations for a Conducting Media 


Equation (V-7) can be solved formalistically as was 


indicated in Chapter III using the Green's function approach. 


That is, 
A(r,t) re ae sah aOR ae 
(V-9) = = fde* fa*r"G (re el] 3", €") P 
Wire) 1 “Pant ' 
, “=—— — (r rt } 
eg, ¢ 


where G(r,t|r',t') is the Green's function and which must 


satisfy the following: 


= 2 Bo LL 
(V-10) {V?-uoo ae ie 


2 eas = ook tt 
aie G(r,t]r',t') = 6(r-r')6(t-t') 
at? 
Equation (V-10) can be solved using Laplace transform 


technigues to yield (see Appendix C) 


& 


~o , or 6. 
a ‘ ee ) Jeet Zery (tt) 2- (Ey? 
(v-11) GF, t]F',t")=7—[e sh sei. enaerceetrneree eae 


r 


Yecety2_ sty 2 
¥ (t-t') (2) 
a ea iy 
oD) Cake te O(t ee =) 
where 15 [x! is the modified Bessel function with argument x and 


u|x| is the unit step function with argument x (U|x|=0 for 
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x < 0 and U[x| = 1 for x > 0). The Pacey should note that 
the above expression was derived assuming that o was con- 

stant with respect to time which, from Table 4, is seen not 
to be the case. Yet, even with this simplifying assumption, 
the author was unable to successfully integrate (V-3) using 
(v-11). This result was included to indicate to the reader 
the difficulties involved in applying the microscopic model 


of the high frequency EMP to conducting media. 
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VI. Final Remarks 

The purpose of this thesis has been to present a model 
for the high altitude, high frequency EMP which could be used 
as a pedagogical tool. As such, the emphasis has been placed 
on developing and outlining the theory rather than obtaining 
state-of-the-art numerical values for threat levels of the 
EMP. However, this author feels that this model has great 
potential for application in this area. The reason for this 
is that this model allows one to see the effect of changing 
a greater number of physical parameters than does the macro~ 
scopic approach which was briefly described in Chapter I. 
For this reason, the following recommendations have been 
made to anyone interested in following up on the work pre- 
sented in this thesis: 

1) An attempt should be made to solve the field 
equations for a conducting media as indicated in 
Chapter V. The approach to this might involve: 

a) finding a numerical solution to the 
integral involving Green's function given 
in Chapter V in equations (V-9) and (V-11) 
or 

b) approaching the problem using a technique 
other than that presented in Chapter V. 

2) If one is able to solve Recommendation #1, then it 
will become necessary to make more precise calcu- 
lations for conductivity than were included in 
Chapter V. One should include the various loss 
mechanisms for the secondary electrons such as 
attachment and recombination. ; 

3) More study should be made regarding the cross- 
sectional line-of-sight area A(R). This could 
include sensitivity studies with respect to 


increasing or decreasing A(R) and the possibility 
of including the n dependence indicated in (III-40). 
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Appendix A 


Gamma Ray Interaction with Matter 


Attenuation Law 

The probability of a photon interacting with the atoms 
of its environment can be expressed in terms of either the 
microscopic cross section or the macroscopic cross section of 
the material which constitutes the environment. The micro- 
scopic cross section (o.) is defined as the effective cross 
sectional area which an atom presents to a gamma photon’ for 
interaction. The macroscopic cross section (yu) is a more 
useful quantity than the microscopic cross section than 
driving a law for the attenuation of gamma radiation. The 
macroscopic cross section is defined as the probability of 
reaction per unit length of photon path and is related to the 
microscopic cross section by the relation 
(Al) u = No, (Ref 1, p 1-61) 
where N is the atom density of the environment. Since oy has 
units of area and N has units of inverse volume, then up must 
have units of inverse length. It is obvious from (Al) that 
pw will depend upon the physical state of the environment at 
any given time, whereas 5 is a fundamental quantity which is 
constant for a given material regardless of the physical state. 


The quantity N can be found from the following relation 
AW 


where N, is Avogadro's number, p is the material density, and 


(A2) N = 


AW is the atomic weight of the matter making up the environment. 


For the purpose of this study an attenuation law with 
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respect to a spherical geometry (see Figure 29) will be 
developed. As shown in Figure 29, the source shall be con- 
"sidered a point source such that photons emanate from it 
egually in all directions. Thus the only variation in the 
photon intensity (photons/unit area) is with respect to the 
distance from the source r. If . represents the number of 
photons which originate from the source and if one assumes 
that there is no removal mechanism for the photons in the 
transport media then the number of photons crossing a unit 
area at a distance r from the source (intensity - $(r)) is 


given by 


$ 
fe) 


(A3) O(r) = 
4nr? 


since 4mr* represents the surface area of a sphere of radius 
x. If a material of cross section KL is placed around the 
source as a shell with an inner radius of r and thickness dr 
(see Figure 29A) then the intensity of unreacted photons 
(virgin intensity) at r + dr will be 

(A4)  ° O(rtdr) = O(r) + do(r) (see Figure 29B) 

In the previous equation dé represents the change in the 
virgin photon intensity across dr because of interaction 
with the material shell. Since » represents the probability 
of an interaction ae unit photon path length, this change 
dé must satisfy the relation 

(A5) dé(r) = -$(r)udr 

That is, the number of photons which interact in the thick- 
ness dr is equal to the total probability of photon interaction 


in the shell (ydr) multiplied by the number of photons "entering" 
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INFINITESTMAL SECTION 
OF THE SHELL PICTURED 
iN PART A TO THE LEFT 
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MATERTAL 


Figure 29. Spherical Geometry for Attenuation Law. 


the shell. The equation (A5) can be rewritten such that 


(a6) SRE} = -uar $ 
If ne pees the gamma seunee as completely surrounded by an 
“infinite number of ‘shells of infinitesimal thicknesses such that 
- the | pHnee most. shell has an inner’ radius of sons and’ the outer- 
most shell has an outer radius of R then the virgin flux at RD 
can be Seance, by integrating the left hand side ar (A6) © 
between ue and 3 (R) and the right hand side between 0 and 


R. That is 


(R) + oR . . 
Sf dao/o = - f pdarf{assume u=p (xr) } 
O,/4mR? oO 
or 
2 'R 
n {2iRi ATR } Se f udr . 
ro) ro) ; 
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or equivalently 


R 
-f ydr 
O(R)47R? _ io 
pgs 5 = oe 
re) 
or solving for $(R) 
R 
® -f ywdr 
(A7) a(R) = ~~ e © 
4nR2 


The last expression above represents the attenuation law 
which is used in this thesis for the virgin intensity of 


gamma photons. . - 


Compton Effect (Ref 16, p 432) 

The primary gamma photon interaction with respect to the 
high altitude EMP is the Compton effect. The analysis of this 
effect is performed by applying the energy and momentum con- 
servation laws to a collision between an electromagnetic 
quantum (gamma photon) and an electron. If EL is the energy 
of the incident gamma, oe the energy of the scattered gamma, 
moc? the rest energy of the electron, and We the total energy 
of the recoil electron (see Figure 6, Chapter IT) then the 


energy conservation law can be expressed as 


BE +mc? = 8 ' + W 
Y ° ¥ e 


Using the substitutions indicated in Figure 6 in Chapter II 
one obtains for the above 


20 ' 
(A8) hv + moc = hv! + We 


Similarly, the momentum conservation laws for this situation 
can be written as 


(A9) hv _ hy! 
aor oe cos 6 + Po Soe B 
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and 


Sa HIN Fite ie 
(A100) 0 = —s7 sin te) Po sin 8 


where the variables used above are those indicated in 

Figure 6. The variable @ can be eliminated by rearranging 
(A9) and (Al0), squaring both sides and adding. The results | 
of this are 


Doe PRU ea) BV Sere 


a2 a2 c2 e 


If one multiplies the above by c* and then rearranges terms, 
the equation becomes 


(All) h? (v-v')? + 2h%vv' (l-cos 6) = Pio? 


Equation {(A8) can be solved for Ws whereupon one: finds 


2 faye ty 2 at 2 oie 3 2 
. {Al12) h* (v-v')* + 2h(v-v Jmo¢ + mvc We 


If one subtracts equation (All) from (Al12) and eliminates We 
and P, through the relativistic relation 

Wo 2 2 pl Bete aa 

e e fe) 

then one obtains 
(A13) 2h(v-v')moc? = 2hvv'(l-cos 6) 
Equation (A13) can then be easily solved for hv' with the 
result that 


het hv 


ioe (1 - cos 6)} 
(@) 


or 
y 


ti 
Hl 


Y {1 + Ev 


or (1 - cos 6)} 
oO 
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Appendix B 


Derivation of E-and B for an Accelerated Electron 


The Wave Equation (Ref 4) 
Maxwell's equations are the starting point for finding 


relations for E and B for a given source. These equations 


are 

(B1) Vx E(r,t) = - ap (r,t) 

(B2) V x B(r,t) = us (x,t) ae Wey Se ire) 4 
eee _ 1” ae 

(B3) V°eE(r,t) = Ey p(t) 

(B4) VeB(r,t) = 0 


where the variables are as indicated in Chapter III. Hence- 
forth it shall be assumed that the conductivity of the media 
surrounding the current source is zero. Because (B4) is 
-always true one can say that 

(B5) B=VxA 

where A is an arbitrary vector at this point. Substituting 
(B5) into (Bl) and assuming that A is such that the curl 
operator and time derivative operator can be interchanged, 


one obtains the following: 


aA (r,t) 4 a, 


(B6) V x {E(r,t) + aE 


Since for any scalar V it is true that VxVV=0, one can make 


the following assertion without loss of generality: 
=, 0A _ 
E + ree VV 


or 
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(By). Se 2 22 ey 


where V is an arbitrary scalar. If one substitutes (B7) 
into (B3) and again assumes that A is such that the spatial 
derivative (V*) can be interchanged with the time derivative 
then the following relation is found: 

p., (Lr t) 


€ 


(B8) Vv (E,t) Hee (Gea (epeh} Ss 
: Vv 


Also, if one substitutes (B7) and (B5) into (B2) and uses 


Pha Wdentitias Ve(VenEV (y mISV oR endl wie Heinen one. 4 


Vvc 

obtains 
ee a ee ee 
(B9) V°A oz bt? bed VAY AtCz aE 


At this point there is still a certain "arbitrariness" with 
respect to A. Since any vector is uniquely defined when both 
its curl and divergence are specified, the arbitrariness in 
A can be removed by specifying. its divergence or the gauge. 
By studying (B8) and (B9) one observes that a convenient 


gauge for this problem is the Lorentz gauge or 


(B10) wAs-+, 2 


Substituting (B10) into (B8) and (B9) one finds the following 


relations for the wave equations of A and V: 


aae 2 — p (r,t) . 
(B11) Ge) ee Si 
c at E., 
and 
2 at 2n iy Ses 
(B12) V?A(r;t) - ty gAte st) = - Hd (r,t) 


The Green's Function (Ref 4) 
Solutions to (B11) and (B12) can be written formal- 


istically using the Green's function approach such that the 
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following is erue:?” 
R(E,t) FE t") 


(B13). = jarr’ fat G(r tle! t!) 


Vi(e,2) 


where G(r,t|r',t') is known as the Green's function and 


which must satisfy the following equation 


2 fae at ' ish ass 
(B14) V2e(x,t]r",t) - 2, VGC tfe' tL geeety) 3 (e-e') 


c ot 
; 92 
One should note that the derivative operators (V? and 5Ez? are 
both with respect to the unprimed coordinates. To solve (B14) 
one can use the Laplace transform technique to remove the time 
derivative from the problem. Thus, if one lets A denote the 


Laplace transform operator with respect to time then the Laplace 


transform of G becomes 


~ foe] 
(B15) G(z,s]rE',t')=A{G}=sfate °* G(z,tl|Flt") 
Oo 
Using the property for Laplace transforms of the pth derivative 


of a function F(t) which is as follows (Ref 8, p 490) 
(p) p Hay ape tan (i) 
A{p‘P? (4)} = sP£(s) - J sP F (+0) 
n=0 


where f(s)=A{F(t)}, one can then write for the Laplace trans- 
form of the second time derivative of G the following: 
2 =e yt ' MAS Sr = pam a 
A{e ASA: eS acae, } = s?G(r,s]r',t') - sG(r,o|r',t') 
_ sG( lett) 


ot t=0 


If one demands that the Green's function and its first time 


19 ohis is just a short hand notation to avoid writing 
both equations. 
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derivative satisfy homogeneous boundary conditions at t=0 then 


the previous equation becomes: 
; : 2 aA, =) ' BE te _ 
(pie) a2 oirek rit’) } = s*e(r,s[r',t') 
dt 


One can find the Laplace transform of the right hand side of 


equation (B14) directly by solving the equation: 


st 


A{6(r-r')6(t-t')} = f date °~ 6(t-t') 
(@) 


which, if one applies the properties of the Dirac delta 


function, becomes 


(B17) MEGS 4). = Sees eo 


Using equations (B15), (B16) and (B17) one can see that if 
the Laplace operator is applied to both sides of equation 
(B14) that this equation becomes 
— oe . ens, =e ae 
(B18) ve G(r,s|r',t') - = Givps|erpe4 = 6(r-r')e a 
c 
where the r subscript to the ve operator is used to emphasize . 


that the derivative is with respect to the unprimed r. By 
translating the origin from 0 to r' and letting R=r-r' one 


can write (B18) as 

= 2 = 2 ' 
(B19) v2 G - = G = 6(R)e 
Cc 


~ 


where the argument of G has been dropped for brevity's sake and 


where it has been assumed that Ve vas Va: (this is easily shown). 
st! 


~ 


' x 
If one multiplies (B19) by et on both sides and lets G=G e 


then this equation can be written as follows: 
x 2 = 
(B20) v2 G- 2G = 6(R) 
R 2 
c 
Since 6(R) represents essentially a point source at R=0, one 


can argue that there is spherical symmetry in the problem and 
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the only spatial dependence which G has is with réspect to the 

aistaues from the source and not on direction. That is, 
= G(IR}) = E® 

Thus for ve G(R) one can write the following 


~ 
~ 


oe aes 2 dG 
VRE oe aR (R dR). 


Using the above and letting G=g/R one can rewrite (B20) as 


2 2 
i d*g _ st 2 = §(R) 
aR? c 
or 
2 2 
(B21) gag 5 Bg = RS IR) 
aR? c? 


If R>O then 6(R)=0 and the right hand side of (B20) is zero; 
. if R=0 then the integral of the right hand side is zero. Thus, 
without loss of generality, one can solve the following instead 


of (B21): 


2 2 
(B22) <i ee eae ee 
aR? -. 


Q 


Equation (B22) has the following general solution: 


e SR/c te et SR/c 


g(R) = Cy 2 


Since it is desired that the fields go to zero as R approaches 


infinity, then C., must be zero in the above. Also, in order 


2 
that the final solutions would reduce to the electrostatic 
solution for V(r,t) it is seen that Cy must equal 77° Thus 
20 


This is: stated without proof. The interested reader 
can carry Cy through the derivation to find V and then match 


the solution to the electrostatic one to find Cy: 
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one has for g(R) that 


1 “SR/c 
g(R) = - qe 
then for G, 
~sR/c 
G = Gf{R) _ =e 
= R 4mR 
and for G, we coe 
aa 
# -st' eae aaa 
G(Z,s{z",t') =e Ge ee 


(B23) 


A general identity concerning Laplace transforms is that 


(B24) 


4n{r-r'| 


A71f{e 8") = 8 (t=) 


-1 ' . , : 
where A denotes the inverse transform. Therefore, if one 


=e | 
lets rat reteor' | and applies equation (B24) to (B23) the 


result is 


G(r,t|r',t') 


G(z,t[r',t") 


(B25) G(r ,t)r se") 


Since the Dirac delta 


A7*iG(=,t EF" ,t')} 


-s (eles hy 


ae 


an |r-r" | 


A~l{e 


} 


i a 
4n|¥-F" | 


Sfe-(et + Jeet hy) 


function is an even function and since 


aaa ’ . alae ' 
e- (er tzce' l= ~fer- (ete hy, equation (B25) can be re- 


written as 


(B26) 


(ef Gage ofl aire one ee 


4n|r-r' | 


Going back to equation (B13) one then has for A(r,t) and 
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V(r,t) the following relation: 


Bie Pe ee ce 
(B27) | | = fe fa°E' sat’ hes cee oe gs 

V(z7t) [r-r' | e P(r rt") 
E and B Solution 


Before one can find a solution for E(r,t) and V(r,t) 

from (B27) one must have an analytical expression for J(r',t') 
and p (et, t'). For an electron with charge e and velocity 
v(t) at the point x(t) the current and charge density are’ 
given by the following: 

J(r',t') = ecB(t)dé{r'-x(t') } 
and 

p(e',t') = ed{r'-x(t')} 
where &(t')=v(t')/c. Thus one can write for the vector and 
- scalar potential the following: 


i 


A(r,t) u_cB(t') 


< 


ae 3s 4 
=q7/4 eUlat 


E-r" 
sfet-(e-EL)] S{E'-(t')} 


eit 


V(r,t) |r-r'| 
v 


Using the properties of the Dirac delta function the integration 


over the spatial coordinates (r') can be immediately performed 


to yield: 
A(xr,t) egntt!) orluceti(t') 
cle eetaut OS ge er 
(B28) V (r,t) R(t") = 
ey 
where R(t') = [r-x(t') | 


One should observe that the only spatial dependence in 
(B28) is in the term R. Further, since the forms of E and B 


given by equation (B5) and (B7) require the use of the spatial 
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operator Me it is convenient to express this operator in terms 
of its effect of R rather than r. Using the Einstein summation 


notation a the vector R can be written as 


R= x,e;5 (i = 1, 2, 3) 


where the Xs represent the coordinates of R in the cartesian 


basis e.- One can also represent the Vi operator by 


or) 
(B29) “ ~ ©4 Ox. 
i 
so that for some function f(R) = F{R(x,)} one can write that 
~  9£(R) : : 
V_f(R) = e5 ax 


= “* df(R) OR 
V_f(R) =e 
xr i aR OXs 
— _ * @R df(R) 
an ae! Gary PRET 
i 
or, using (B29) 
= _ op d£(R) 
(B30) Vf (R) = VR sR 
: 7 1/2 . — 
Since R=(x;x,) one can write for VR that 
VR # e, @= (x. %,)1/7 = ey 272545 
} Ox. ivi j 172 
2(x.x.) 
. Loi 
or 
x.e. = A 
Rr- J 3 = ae 
(B31) VR V3 = =n 
(x.x.) 
Li 


aw 
where, because of the definition of R, n is a unit vector 


21yy this notation a repeated index signifies a sum. 


Thus a,b, = . a,b; = ajb, + anb, te ae atectrdes 
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pointing from the source to the observer. Thus one can write, 


using (B30) and (B31), that 


. A 
(B32) Lee *~ Nn oR 


Using the relation for E given in (B7) as applied to (B28) 


and letting Ve be as given in (B32) one finds the following 


for E(r,t) . 
Biro) 2 So yee (eee oie ee 
4 eR (t') c 7 
bys R(t") 

a(t! + S24 -t) 
Rit ) c . 

€ rire pl ortitae aR 
B(t')uc ' ' _ 

est ets See Vv {é! (t! eRe) ) -t) a(t tat )/c t }) 
R(t") 


where 6'(x) denotes the first derivative of the Dirac delta 
function with respect to its argument. After terms are re- 
arranged and the indicated derivatives performed the above 


expression for E(r,t) can be written as 


“nw 


(B33) E(¥,t) = 7=— sat ences é{t' + Be -t} 
ame. R2 (t') c 
+ RENE (B(t') =n) Go “nd ~t}) 


where the fact that Hoe, = ty has been used. In order to 
evaluate the integral in (B33) it will be divided into the 


following two parts: 


fat! (a 
Re") 


(B34) I 


1, R(t") _ 


(B35) z 


1 = a ' RCtT JS 
2 = fat" e7ENS (B(t')-n) {S'(t' + ns t)} 
It should be noted that the argument of the delta function 


in I, (and.T,) contains a function of the variable of integration 


100 
“43 


GNE/PH/74-7 


rather than simply the variable itself. In order to evaluate 
this particular situation the following general integral will 
be considered: 


I= fg(t') 6{£(t') -c}dt' 


where g and f are arbitrary functions of x each of which has 
a dagepae*- and first derivative and where c is some constant 
with respect to the integration variable. If one lets £(t")=u 
and then makes a change of variable in the integral such that 
t'=£ )(u) then " 2 


_ at F ee > One, 
du = ae dt' + dt' = df/ate 


and I becomes 


_ , att?) 7 
I = t df/at® 6 (u c)du 


or using the delta function properties 


nas Ste) 
adf/dt' lu=c . 
or 
(B36) I= fg(x)é6{f£(t')-c}dt' = 2 g{t,") 


i df/dt'|£(t')=e 
where the t,' are those values of t' which satisfy f(t')=c. 
Applying (B36) to (B34) one can see first that £(t')=t'4or(t") 


which implies the following: 


di _ 1 aR 
(B37) gev = 1 + Caer 
If one represents R by (R-R) 1/2 then 4 Soy can be written as: 
LGR _1 a gigyl/2_1 Rak 
c dt' c dt’ c R dt 


or using the definition of R and letting nk 


22nor example if f(x)=u, then x=f 1(u) where f ! is the 
inverse function of f. . 
Toi 
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. 1 dR 
(B38) = ar 


Thus one has for .(B37) 


at. 


date =] -n-8(t')=Ek 


(B39) 


Using (B38) and the relation defined in (B36) one finds that 


(B34) can be written as follows 


n 


7 _ n 
TO re ee ee 1 
Po er?(et) Jer + BED et 
or 
- n = i 
(B40) qT ~ BP peng a R(t') n 
KR*(t') jtt=t - ——— Pe: asa 
c eR? 


The brackets with the subscript "ret" means that the quantity 
which they enclose is to be evaluated at the retarded time, 


ete 8, 
Cc 


The expression (B35) and I, can be evaluated similarly to 


2 


if one makes the following substitutions: 


1 
f(t") =e + RED 2 ye th = eu) 
(B41) Gers aEvaET = oo ( as before) 
oa! “~y = io 
g(u) = se oes ae 
R(f ?) 


Then one has for I, that 


=i 5 SO). 6 (y-t)au 


I, - 6! at/de 


The above can be integrated by parts on the derivative of the 


delta function to yield: 


Poe -+ PRGRY 


: = 
2 at' K du | 


This can then be integrated in the same manner as I. to yield 


ali 
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i 


dt'd su) }| 
du dt' K u=t 


ae 
> ae c ( 


or uSing the substitutions indicated in (B40) 


=. |L. d_ /8=n 
sate) > aa Ke dt! CR ee 


2 


The equation (B33) for E(r,t) can be written in a closed 


form if one employs (B40) and (B41). The results of this are 
— e n see | n- B 
ame, cR2 Ke dt KR eee 
To perform the differentiation - * of one can use the following 
identity presented here without proof or 4ustification?* 
Lar, Pe oe 
(B43) cs att 7 R 1 (n-8)n-B} 


Using the substitution indicated by (B43) and performing the 


necessary algebra, one finds the following relation for 


E(r,t): 

—_,— e n n dst 8 Ll B 

E — ait els pet ae fn Se rn ley Pantha 
(Beth: EAE) 4mé | 2p2 cK act (eR) 2,2 GK GE" Ce sit 


If one had performed similar calculations to those performed 
up to this point on B(r,t) where B=VxA then the following 


relation would have been found for B(r,t): 


ne ecu aD R * 
(B45) -B(E,t) = pe (2 + & SHA) «x 0] 
K?R? 


One can easily show using (B44) and vector algebra that the 


following relation exist between E and B: 


(B46) B(E,t) = 2 fy ok Beet 


2376 verify this relation one should realize that 
r=x(t')+R(t') and that for a fixed observer 2£=0. Also use 


the definition of n=R(R*°R) me 
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In order to obtain a more explicit form for E(r,t) 


than that given by (B44) the following relation is noted: 


lad ee ge 

Caer SR) Se ae a " B)R| 
(B47) 1 foie aie an 1 dR 

Cc Gor (KR) = Bins B) a8 at? + c GR (a-n ne 8) 
where 


(B48) B= SE 


Equation (B47) can be rewritten using the substitutions 
indicated by (B38) and (B43) such that 

ld 
c 


d- (eR) = 8? - (Ben) - B(n-B) 


(B49) 5 


Therefore, if one uses the substitutions given by equations 
(B48) and (B49) and rearranges the terms using vector algebra 


then equation (B44) can be written as follows: 


ae ee n-B) (1-|B}2) 1 : 
re eee Te. { ie base a [Zale B) xB} oa! 
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Appendix -C 


Derivation of the Green's Function for Conducting Media 


If one assumes that a conducting media obeys Ohm's law 
and that the conductivity is uniform throughout space and 
constant with respect to time then the following wave equation 
for the vector and scalar potential can be derived (see 
Chaptet V) 

5 3. 22 A(r,t) —U5 (x,t) 


(cl) {V?-po = -= 
vi Bt og? at? 


V(x,t) | tL p(t) 
Vv 


The Green's function associated with this must satisfy the 


following equation 


8. 2 8) gE, e]F,e") = 6 (E-F') 8 (t-t") 


(C2) {V?-u_o —— 
ce? ot? 


vot 
One can solve (C2) in exactly the same manner as indicated in 
the second section of Appendix B using Laplace transform tech- 
niques. Doing this one finds the following for the Laplace 


transform of G(r,t|r',t'): 


(C3) G(Z,s|E',t!') = — e e 

It will be necessary to rearrange some of the terms and use 
some transform identities in order to get (C3) into a form 
which can be found in a table of Laplace transforms. The 


following identity shall be employed (Ref 8, p 491): 
(C4) A +{£(s)e } = F(t-t') 


where 


A7?{£(s)} = F(t) and F(t) = 0 for t<0 
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Therefore it is only necessary to find the inverse transform 


of the following for (C3): 


Lr 
es oy, s (_— + s} 
(C5) g(s) =e Vv 
If one lets K=¢— then one can write that 
Vv 
s(s + ae s(s + K) 
fy 
or 
s(s + K) = s* + sk + (5)? - (5)? 
or . 
(C6) s(s + K) = (s + £)?-(§)? 


Using (C6) one can then write the following for (C5): 
xr / K K 
= (s + ot = oa 
(C7) g(s) =e 
“If F(t) = A”'{h(s)} then the following identity is true 
(Ref 8, p 491): 
(c8) .  AMn(s - b)} = e&* £(t) 
where b is a constant with respect to the transform. Thus, 
if one lets h(s) be the following: 
x 
Gv st = (5)* 
h(s) =e 
then from (C7) 
Xr K K 
~@ 7 (s 4° - G)? 
h(s-b) = g(s) =e 


where, in this case, b=-5. Thus, using (C8) one finds 


(c9) A *{g(s)} = A ‘{h(s-b)} =e 2 


The last term in (C9) can be written as follows by adding 


0 to both sides: 
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~e } 


The inverse transforms of each term on the right hand side of 


(C10) are common in tables of Laplace transforms. 


a table the following is found: 


= s?~ (8)? 
cv 2 Ke jk /+2-(F)? 
-1 2c 1!'2 c 
(C11) A fe } = 
ry 2 

leas 
r 

+ § (t-2) 


Using such 


Olt 


where 1, |x| is the modified Bessel's function with argument 


x and U|x| is the unit step function with argument x. 


stituting (C11) into (C9) one finds that: 


=e t 
2€ or fe} fi2 (ty 2 
24 ¥ 2e,c Tie. . te) | 
(C12) A i{g(s)} =e { ——__________________— 


/-(2)? 


+ &(t-5)} 


Sub- 


ult-=| 
Cc 


Finally if one applies (C12) with the identity given in (C4) 


to (C3) the result is as follows: 


aH} 


G(z,t{z',t') = AT {G(E,slE',t")} 


10) 
“Ze, (tt!) g 
ee , wu 1,|2e, Wet") 2-(2)?| 
(C13) G(r, t[F',t") = ple (ge meee een emer 


¥(t-t')?=(E)? 


x Ult-t'-E| + 6(t-t'-=))} 
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